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1,0 IKTRODUCTION AND SUMMARY 


The Blade Loss Transient Dynamic Analysis Program comprises three tech™ 
nical tasks, namely: 

0 Task I ” Survey and perspective of problem areas in engine transient 
dynamics and selection of specific problem areas for analysis 

0 Task II “ Analysis and development of analytical techniques of 
engine transient response from dynamic stimuli 

• Task III ™ Experimental verification of the analytical method by 
correlation of analysis results with test data. 

The ultimate objective is the development of a computer code for the 
analysis and prediction of turbine engine transient response from the loss of 
a blade. This has been necessitated by the growing awareness that engines 
must be lighter, more rugged, and rim at higher speeds to achieve maximum 
efficiency and economy. The implications of blade loss, and subsequent 
dynamic loads, on engine integrity and performance become obvious. In order 
to tolerate unbalance loads, engine response must be predicted accurately so 
that an optimum structure can be designed with adequate bearings and dampers. 
Thus, damage will be minimized, in the event of a blade loss, without unduly 
penalizing the engine with regards to weight and fuel economy during normal 
operation. 

From the results of Task I, a survey of field experience of blade loss 
and its consequences in turbine engines, critical problem areas were identi™ 
fied. The most critical problem areas with respect to probability of occur- 
rence and subsequent damage to the engine were selected for analytical con- 
sideration in Task II. 

The critical problem areas may be classified by fundamental mechanisms 
of resulting damage: 

0 Deflection Dominated - This is primarily a rub between rotor and 
casing or rotor-to-rotor 

0 High Loads and Stresses - Such as vibratory and transient loads. 

The analytical approach chosen to develop a transient dynamic analysis 
computer program in Task II is the component element method. This method is 
essentially based on modal synthesis combined with a central, finite-differ- 
ence, numerical-integration scheme. The methodology leads to a modular or 
building-block technique that is amenable to computer programming. It also 
presents a simple way to model engines of various structural configurations. 



In order to verify the analytical method, Turbine Engine Transient 
Response Analysis (TETRA) developed in Task II, it was applied to two blade-- 
out test vehicles that had been previously instrumented and tested. Compari- 
son of the time-dependent test data with those predicted by TETRA is the sub- 
ject of Task III. Finally, these comparisons led to recommendations for 
refinement or extension of the analytical method to improve its accuracy and 
overcome its shortcomings. 

The subject contract has been completed, and the results are documented 
in the following reports: 

• Volume I contains the results of Task I. 

e Volume II contains the theoretical basis for the method of analysis 
of Task II as well as results of Task III. Included are the tech- 
nical details of Task II (the application of the component element 
method to the development of TETRA). This includes the development 
of working equations, their discretisation, numerical solution 
scheme, the modular concept of engine modeling, the program logical 
structure, and some illustrated results. It also contains the ana- 
lytical predictions for two blade-out tests as well as the experi- 
mental verification required in Task III, Described are the blade- 
loss test vehicles (rig and full engine), the type of measured data, 
and the engine structural models. It also presents a comparison of 
the predicted response of these vehicles with sudden unbalance as 
calculated by TETRA with the appropriate measured dynamic data. 
Finally, the shortcomings, inadequacies, and recommendations for 
refinement or increased capability of program TETRA are reviewed. 

• Volume III is a User’s Manual for TETRA, It contains program logic, 
flow charts, error messages, input sheets, modeling instructions, 
option descriptions, input variable descriptions, and demonstration 
problems. The process of obtaining a NASTRAN 17.5 generated modal 
input file for TETRA is also described with a worked sample. 

Each volume is divided into sections treating major topics in detail. 

In Volume II, for instance, Section 1.0 treats the fundamentals of the compo- 
nent element method in a clear and unified fashion and shows how physical 
elements and modal components combine to produce a representation of the whole 
engine. Section 2.Q of the same volume is much more detailed due to the many 
small items that must be coded in a computer program. Here the mathematical 
representations of the many modules are developed in detail: modal mass, 

modal stiffness, spring, dashpot, rub element, gryroscopic forces, unbalance- 
force excitation, time-integration scheme, total displacement, etc. Because 
of the modal approach, assembly of the engine is seen in the way one makes the 
input sheets to model the entire structure. 

The volumes are written so that each is self-contained; therefore, some 
repetition is unavoidable. 
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2.0 MODAL SYNTHESIS AND THE LAGRANGIAN IN THE DEVELOPMENT 
OF THE COMPONENT ELEMENT METHOD 


2.1 • SUMMARY 

The fundamental principles of the component element method are explained 
and illustrated. The component element method is based on the modal synthesis 
approach; the response of a complex structure is constructed in terms of the 
uncoupled normal modes of the main structural subsystems with appropriate 
connection forces introduced by physical connecting elements such as dashpots 
and springs. Lagrangian dynamics is introduced as a simple way to obtain the 
equations of the combined system. This approach also allows derivation of 
additional connecting forces such as those from squeeze films, rubs, and 
others in a general way. 

The building block or modular approach follows directly from the logical 
structure of the modal synthesis approach. By representing the rotor motion 
with its nonrotating normal modes and the additional displacements and veloc" 
ities due to rotation, it is found that the g 3 nroscopic forces are reproduced 
and can be treated as "connection" forces tdiich couple the motions in two 
orthogonal planes. 

It is also shown that the classical rotor equations neglect so-called 
parametric excitation introduced by unbalance. These parametric terms are 
contained in the harmonic or time-dependent coefficients of rotor angular dis- 
placement, velocity, and acceleration at locations t^ere there is an unbal- 
ance; they are comparable to hatmonic variations in rotor moments of inertia. 
Omission of these terms is found to be justified since they are proportional 
to the square of the unbalance eccentric radius. 

This section mainly clarifies the component element method so that its 
fundamentals are more easily understood. The actual working equations used 
in the development of the TETRA program are given in Volume III. 


2.2 BACKGROUND AND INTRODUCTION 


Carpenters build houses with 2 x 4's, rafters, concrete slab, and nails; 
we make turbine engines with blades, disks, combustors, rotors, and casings. 
Therefore, it is necessary to develop a method for predicting the dynamic 
stresses and vibrations of a complex structure in terms of the behavior of 
individual component parts. Modal synthesis is such a method of calculating 
structural dynamic response. Under various appellations such as transfer 
function (Reference l), impedance (Reference 2), admittance (Reference 2), 
component mode (References 2, 3, and 4), or component element (Reference 5), 
which are esoteric in different technical disciplines, modal synthesis is 
an elegant but simple concept that may have been developed as far back as 
the 19th century . Mathematicians associate it with the solution of partial 
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differential equations by expansion in orthogonal functions, such as Fourier 
Series, and call it Rayleigh«-Ritz (Reference 6), or Galerkin*s Method when 
the functions are not orthogonal (Reference 6). 

The principle is that the dynamic response of a complex structure is con** 
strutted by the ordered concatenation of the characteristic responses of its 
parts. These characteristic responses are also called the normal modes of 
that component and are calculated independent of the influence of other com* 
ponents. A component may be one simple beam, or it may be a conglomeration of 
several elements. An example of the latter is a compressor rotor which is 
made up of rotor, disks, blades and seals; the normal mode would be that of 
the complete compressor rotor behaving as a single entity. 

These various components are assembled to form the complete engine with 
connecting elements at bearings and other support points. These connections 
not only couple the components but also modify their boundary conditions. 
Because the modal characteristics follow the same form and connecting elements 
are expressed by their damping and spring characteristics, the basis for a 
modular representation of complex structures is a natural consequence. 

To obtain the proper and ordered concatenation of the components modes 
through the connecting elements, as well as to assure the correctness of the 
dynamic equations, Lagrangian dynamics is employed. By starting from the 
fundamental concepts of energy and Euler-Lagrange variational principle, the 
equations of each component and the proper connecting terms can be derived 
without ambiguity. 

The resulting equations of the assembly can be solved by several tech*” 
niques such as Runge”*Kutta or the "difference” methods for transient response 
calculations. The marriage of the modular representation of the modal syn*“ 
thesis method with the central-difference, numerical-integration solution is 
the essence of the component element method. 

This report will detail the fundamentals of modal synthesis and the 
Lagrangian approach for the analysis of structural dynamic response in gen- 
eral and turbine engine dynamic response in particular. 


2.3 MODAL SYNTHESIS 


It is known that any function which satisfies a linear partial differen- 
tial equation (PDE) and the appropriate boundary conditions is a solution, and 
the linear combination of such functions is also a solution. These functions 
have been obtained by intuition, trial and error, and (more directly) from 
the solution of the homogeneous PDE (e.g., with the right-hand side equal to 
zero). The latter is called the Eigen value or characteristic-value problem. 
In dynamics and vibrations, the Eigen values or characteristic roots are 
called the natural frequencies or critical speeds; the deflection shape corre- 
sponding to each natural frequency is the mode shape or Eigen vector. The 
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solution of the full PDE from a combination of these modes is called the 
method of superposition of normal modes or modal expansion- 

Modal synthesis carries modal expansion one step further by utilizing 
the normal modes to calculate the response of the structure when additional 
constraints are imposed- 

These mode shapes are also called normal or orthogonal modes because the 
product of any two mode shapes when multiplied by the mass and integrated over 
the volume is zero when the modes are not identical and not zero when they are. 
This is the principle of orthogonality of normal modes. This principle is 
especially interesting because a purely mathematical concept leads to the 
principle of the conservation of energy. 

The modal synthesis with a simply supported beam under an impulse load 
and with a spring at z = z* is illustrated. The PDE and appropriate boundary 
conditions are 


EIu^''^(z,t) + uu(z,t) ® “ Ku{z*,t) + FQ(z,t)5(z - Zo)6(t) 
u(o,t ) « u(L,t ) » u**(o,t) « u”(L,t ) * 0 

where the Dirac delta functions describe the force to act at a point Zq and at 
a time zero. 

The homogeneous equation is obtained with the right-hand side equal to 
zero, and the well known Eigen values or frequencies and mode shapes are: 

0 )^ ® nir^^EI/uL^ 

UqCz) “ sin nirz/L 

Now, assume the solution to be: 

u(z,t) H ^ ui(z)qi(t) 
i 


q^Ct) * unknown participation of i-th mode or generalized 
coordinate 

Substituting these into the full PDE and simplifying gives: 


El 




+ yq 


sin 


nirz 


-K V q &in + F ( z , t ) 

^ L o o 


O O 
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Multiply this with = sin 6q^ and integrate over the span, then 

obtain the partial derivative with respect to dq^. Thus, 


*^n * ‘^^’n f ^ iSiS. dz » -K ^si 


ir if 

mrz . mirz ^ 
in ~ — sin q + 

L L ^n 


* F U .t) ,in^ 

o o L 


By using the principle of orthogonality, the summation sign in the left- 
hand side is dropped. However, the summation in the right-hand side does not 
vanish. Thus, 




. 2 
sin 


n-irz 

L 


dz » -K 




sin 


nffz 

L 


sin 


mirz 


q F 
^n o 


(z ,t) sin 
o 


mirZQ 

«_ 


The first integral is the n-th modal stiffness, and the second integral is 
the n-th modal or generalized masSe In the more standard notation, we have the 
set of ordinary differential equations (ODE) with ”t*' the only independent 
variable: 


M q 
®n ^ 


K q 
®n ^ 


K q 
mn^n 


F 

o 


sin 


mirzo 

L 


M 




uu u 
mm 


dz 



. rairz , 
sin dz 



0,i = j 
=» 0,i j 


where: 


K =» 
% 



mtrz , 

sin 



M 

mn 


K 

®n 


K sin 


if 

mifZ 

L 


sin 


ntrx* 

L 
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It can be seen that these modal equations have coupling terms through 
the support spring. This spring actually alters the whole beam. The normal 
modes of this new beam with the spring at 2 » z* can be obtained in terms of 
the modes of the original beam by solving the homogeneous equation of the new 
beam equation. 

This example is a simple illustration of modal synthesis. Now suppose 
we have two simply supported beams that are parallel to each other and con- 
nected by a spring. This is illustrated below: 



u(0,t) » u(L,t) » u»(0,t) “ u*’ 
(L,t) “ 0 

v(0,t) ® v(L,t) « u”(0,t) * u** 
(L,t) « 0 


The governing equation of each beam is similar to that of the single beam. 

The difference is that the spring force at the connection points is propor- 
tional to the relative deflection and, hence, mode shapes of both beams. The 
PDE of both beams are: 


* Fo(zo,t) “ K[u(z*,t) - v(z*,t)] 

ElyV^^ + UvV « K[u( 2 *,t) “ v(z*,t)] 

Now writing the beam displacements in terms of their simply supported 
normal modes 

u(z,t) ~ 
v(z,t) 


Zxii(z) qi(t) 
Svi(z) ?i(t) 
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and multiplying the first by (virtual modal displacement of the first 

beam) , the second by (virtual modal displacement of second beam) , then 

integrating over z we get: 


«<lk|qnMkn «ln^ 2C<“1*<U * ••• «n «ln) “ + ... v* = FoUit5qij 

5^k l^n^kn ■*“ " *"k 2[(“1*<?1 + • • • % qn> “ (v* ?! + . • • v*P„)]| = 0 


Taking the derivative of the first with respect to Sq^ and the second 
with respect to 6Pn, we have the ODE of the connected beams; thus: 

^Omn + ^<?n + ^il^l ^n2<l2 “ ^1^1 " W?2 ••• ^?ln 

^^^m^ ■*' ■*■ ^1^1 ■*■ ^2^2 ••• “ ^1*11 ~ ^2*12 




where 


k « oj M M 

°n 


Li 

/ u U U 

u ra n 


dz 0} ® (rtr) -mi j~ 

W tA 




K =* M 
™n “n “n 


M ® i u V V dz 
m n 


0 ) ® (rir) 


•V5 


^’Vi’^’^n’“™n’“®n 



0 , m«n 

® 0, mfn 


= ki%(z ) u^i(z ) 

=• Kv„*Un*J k = Ku^*v„* 


s 



One can see that the modal equations of each beam have coupling terms in*™ 
volving other modes of that same beam and also modes of the other beam. 

The last example illustrates the modal synthesis approach using the nor^“ 
mal modes of the component beams and the characteristic of the spring connec*™ 
tion. From this illustration, the analysis can be easily extended to more 
complex structures. 

However, in application to structures with rotating components as well 
as static elements, formation of the generalized equations is not clear and 
sometimes ambiguous. To avoid this difficulty, the concept of the Lagrangian 
will be described. 


2.4 LAGRANGIAN DYNAMICS 

It is called Lagrange’s equation in the U.S.; however, it is known as 
Euler-Lagrange in Western Europe and in Moscow as the Euler-Ostrogrodski 
equation. For consolidation, it will be referred to as the Euler-Lagrange- 
Ostrogrodski (ELO) variational equation in this report. 

The ELO is a relationship that must be satisfied in order to minimize the 
kinetic and potential energies of a system. It is written thus: 

- iii - F 

dt 35 35 5 

where: 

L * T *“ V, Lagrangian; T - Kinetic energy, V = potential energy 

^(t) ® Dependent variable or generalized coordinate 

F^ Generalized force in ?(t) 

t ® Time; % ® d^/dt 

This is another statement of Newton’s Law that force equals mass times 
acceleration. It is a more general relation since it allows the force rela™* 
tions to be obtained in nonphysical space such as the space of the modal or 
generalized coordinates we have mentioned earlier. 

In the section on modal synthesis, we began with the equation of motion 
of the structure, the PDE, and proceeded to obtain solutions. By contrast, 
the ELO begins with displacements, proceeds to obtain the kinetic and ppten* 
tial energies, and finally develops the equations of motion or PDE. 

One can approach the Lagrangian using the physical displacements which 
would result in the PDE with the spatial coordinates and time as independent 
variables. This is the standard procedure for obtaining the PDE of fluids or 
elasticity such as those in plates and shells and the simple beam discussed 
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earlier. Or the Lagrangian can be expressed directly in terms of generalized 
modal coordinates (e.g,, write the physical displacements and velocities in 
terms of the normal modes). The latter is the typical approach to analysis of 
complex structures or problems which involve the interaction of a structure 
with a fluid, as in aeroelasticity . 

A simplification can also be made on the ELO operation on kinetic energy. 
This is particularly helpful for an articuiacdi j-^.-^care with a complex 
motion and for rotating structures such as rotors or propeller blades. Recall 
that the kinetic energy is the sum of the squares of the physical velocity 
components in the three cartesian axes and that these physical displacements 
may be expressed in terms of modal coordinates or other coordinates fixed to 
the structure. 

It can be shown that the physical displacements measured from a fixed 
reference frame can be written as: 

X - x[c(t)] 
y = y[5(t)] 
z = z[c(t)] 


where 5 assumes all the generalized coordinates and x ( 5 ), say, means a func~ 
tion of all the generalized coordinates. Then, the kinetic energy for a body 
is : 

T = 1/2 /(i2 + y2 + g2) 

The ELO operation on the kinetic energy can be shown to reduce to: 


d 3 T _ 3 T _ 

dt 35 “ H 


3x 

3C 


3C 


z —) 

35 ^ 


dm 


When the final equations of motion will be linearized, the above relation 
will allow separate linearization of the accelerations and directional deriv- 
atives. This will drastically reduce algebraic manipulation and the chances 
of errors in reaching the final result. 

The potential energy of connecting elements can be formulated for each 
configuration in terms of displacements at its ends: A, B. For example, a 

simple linear spring has the energy: 
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The modal energy is much simpler. From the conservation of energy, the modal 
potential energy is simply: 


1 2 
T“(ij) 


M. . 

2.J 


^ij5 w(ij) 



i«j 

in 


In the case of rotating structures such as a turbine rotor, the normal 
modes can be obtained with or without the effects of rotation. With the 
latter, the effects of rotation can be added separately by constructing the 
equations of the rotating rotor with the static normal modes in two perpendic- 
ular planes and the additional displacements due to rotation. Following the 
Lagrangian approach discussed above, the physical displacement of a point on 
the rotor can be written down by considering the geometry of motion. The ELO 
.operation on the kinetic energy will yield the inertial equations of motion 
with the additional gyroscopic terms. Neither modal potential energy nor the 
potential energy of the connecting elements are affected by the rotation. 


2.5 DAMPING AND OTHER EXTERNAL FORCES 

So far, we have shown derivations of the equations of a conservative 
system without external forces. Damping and external forces or even other 
springs can be derived in a similar manner from physical forces with the prin- 
ciple of virtual work. The principle of virtual work may also be considered 
a special variational principal like the ELO. 


2.5.1 Damping 

In the case of viscous damping, Rayleigh *s dissipation function can be 
constructed. It is normally written: 


where : 

D Dissipation function 
C ^ Viscous damping coefficient 

z « ~ , Physical relative velocity between ends of the damper 
or dashpot. 
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If Che dashpot connecCs Cwo structures, a and b, whose norrnal modes are 
known, one can write the dissipation function, say for damping in Che 2 direc”* 
tion: 

2 2 

D « “ C (2 (x ) *“ 4 C (z “ Zv ) 

z a D z a 0 

and Che modal force is: 


F 





• *v 3 f' * 
- Z5 ) — r (Zg 

3C 


. * 

Zb 


) 


where 5 takes on all the possible generalized velocities. The terms Z3 and 2^3 
are expressed as the sum of the normal modes of structures a and b whose modal 
participations are denoted by c(t). 


A separate dissipation function is, of course^ defined for each dashpot. 
These damping forces introduce coupling between modes as well as between com- 
ponents. 


Distributed damping may also be constructed, although this is not always 
practical. However, modal damping has been used frequently. It is an espe- 
cially distributed damping proportional to either the mass or the stiffness 
distribution or both. This proportionality to either of these Cwo distribu- 
tions allows one to take advantage of the orthogonality of normal modes. 
Therefore, the modal dissipation function may be written simply as: 


1 . . . , 

*2 5 ^ij^i^j Mass-proportronal dissipation function 


2 « 


Sci£fness~proporCional dlssipacion function 


and Che damping forces are: 



a 

3qi 

C M, .q . 

3Pk _ 

n K Q 

3qi 

Q fs., . q . 

J 


noting that Mij and Kij are nonzero for ij only. 
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2*5.2 External Forces 


Physical external forces which are functions of location and time act 
along directions with respect to the fixed reference frame. To obtain the 
generalized forces which act on the generalized coordinate equations of mo" 
tion, the principle of virtual work is employed. This principle was used 
earlier in the exposition of the concept of the modal synthesis. The virtual 
work of a physical force is the product of that force times the virtual dis- 
placement along the action line of that force. Consider a force Fg parallel 
to the z axis. The virtual displacement is 6z, and the virtual 6W work is;. 

5 W « - F26Z 

The generalized force in the generalized coordinate is then: 

F - = - F 

c TuU 2 TuT) 35 

A physical spring or dashpot may be also considered as an external force. 
Applying the virtual work principle will yield a generalized spring force 
identical to that obtained in the ELO operation on the spring potential 
energy. 

Recalling the Lagrangian and the ELO and combining these with the dis- 
sipation function and the virtual work principle, the complete equations of 
motion of any complex structure may be obtained with a variational principle. 
Thus: 


F ^ 8L _ 3L ^ 3D „ 3UW) 

5 . dt 35 “35 35 ~ TTSTT 

where: L, D, and contain all the total energies, dissipation functions, and 

virtual work in the whole system, and E assumes all the generalized 
coordinates. 

The resulting fully coupled equations define a matrix whose order is the 
total number of all the generalized coordinates of all the component struc- 
tures, and the only independent variable is time. 

We have shown a method for deriving the equations of motion of single 
structures and more complex structures made up of many components. From these 
governing equations, one can establish the ordered concatenation of the indi- 
vidual elements and the proper treatment of the interconnections. Each mode 
represented by modal mass, stiffness, and mode shape is a module, and a single 
component which is a physical module can be represented by its modal modules. 
With the modal building blocks and the modules of connecting elements and 
external forces, the model of a complex structure can be constructed. The 
solution is the response of the entire system given in terms of the response 
of the normal modes of its components. 
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2e6 APPLICATION TO TURBINE ENGINE DYNAMICS 


A turbine consisting of a casing and a rotor is sketched in Figures 1 
and 2. The casing is connected to ground by springs and dashpots; the rotor 
is supported by the casing with similar connecting elements. 


Because the casing 
is simply: 


is not rotating, its Lagrangian (without connections) 


L 

c 



V « 
c 


1 • • 
-r M, .q.q . 
2 


c 

h 



c 

h 





» « 

M. .P.P. 


c 

V 


tu. . ..M. .P.P , 
(xj) ij J 


c 

V 


where the subscripts 


h Denote horizontal plane 
V Denote vertical plane 


and 


q£ ® Generalized coordinate in horizontal plane 
Pi « Generalized coordinate in vertical plane 


Mij ~ Modal mass; 


t 0, i=j 


0, i?^j 


tuij ® Modal frequency 

The potential and dissipation function of all the spring supports and 
dashpots is simply: 




springs 


springs 




c r 


*1 

z__ ] ^ 




dash 


dash 
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ioirntry 




o: 


[•y»} 


;oi 



where: 


'T*^hi(x*)qj^ yc* “ y^vi(x*)Pj 




%i. 

r 



Since we are using nonrotating modes for the rotor, which have no gyro- 
scopic effect, we will first write down Che displacement components of a point 
on the rotor which is r from the centerline* These consist of the stationary 
mode displacements and those induced by the rotor speed as. well as the small 
rotations from bending slope. Let: 


z(x,t) 

x(x,t) 

t • (at ff) 
a 


Total displacement with respect to ground-horizontal plane 

Total displacement with respect to ground-vertical plane 

Total displacement with respect to ground induced by 
bending, in both planes 

Angular displacement about axis 

Rotor speed, assumed constant 

Stationary modes - horizontal plane 

Stationary modes - verical plane 


a ^ 3z/33C « 3Z/3X ; 6 “ 3z/3x « 3z/3x 


We have not shown Che subscript, for brevity* Also, for conciseness, we 
will express the stationary modes as z and y; a and 0 are the x derivatives. 
It should be understood that these are made up of sums of the normal modes 
h^cji and ^i^i* 

From the figure of the geometry of motion given below, we can write down 
the total displacement components with this sequence of rotations: 

. 1* Rotor rotation: a 

2. Slope rotation about the horizontal fixed axis: a, o 

3, Slope rotation about the vertical fixed axis: 0, 0 
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thus : 


z = 2 + r cos ® cos 

y = y 4- r sin Cl sin 9 cos + r cos 9 sin ^ 

X ® sin ® cos Q cos ^ + r sin ^ sin 


The modal potential energy is not affected by the rotor revolution. So, 
the rotor Lagrangian (exclusive of connecting elements which were already con- 
sidered with the casing) is: 


L 

r 



^ 2 
y + 




r 


w.. M. .q.q. 
(ij) ij 1 J 



2 


2 


(ij) 


M. .P. 





To obtain the kinetic or inertial terms, we operate with the ELO on the rotor 
kinetic energy: 


f 


6 


dt ae 



r 


dm 


with 0 =* 2, y, e. This is done so the results can be compared with those 
normally derived for rotor critical speed. And to obtain the equation in the 
modal generalized coordinates (5 ® <lii^*i)j we proceed as follows 

f 11 

dt 35 W dt ■5F" j dT ZmJ dC 

S' ' B 


and: 


2 * E hi qi 



i 


y » £ vi?i 


3vi 

* ■ r 'i 
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Proceeding in this manner, the kinetic or inertial rotor equations in 
the S * X, y, a, 6 coordinates and then the modal generalized coordinates 
^i> ^i found to be: 


I.- Phvsical Disnlacements 


J* I z a r cos ^ I dm 
fy =* I y “■ a r sin ^ | dm 

f\9T^si 


1 2 
» <1 

2 


2 

r i 2 2 • 2 2 * 22 ) 

ffl ® / ler sin i|> + 290r sin cos - Sr sin ^ cos ^ + 2aCr sin > dm 


.. 2 
er 5 

.. 2 . 


A ^ *22) 

ir cos > • 2aar‘'sin ij; cos - er'^sin Tj; cos ^ - 2anr cos \|i f dm 


In Modal Generalized Coordinates 


f* \ / 2 1 f ,2 I I • 2 ^ t I •• 2 ^ 

fqi « / jCh^hj ^ r h£hj cos t|i) qj - 2h£hjqjar sin ^ cos hiVjPjr s.in ^ cos tp 

t ,•» 2 2 2 , 

2h£VjPjQr cos “ a h^r cos dm 

r i 2 , f 2 ft* 2 ^ 11^ 2 ^ 

fpi * J + r V£Vj sin Pj 2v£Vj PjQr sin^costj; - Vjthjqjr simjicosii^ 


2 2 2 

2v!h!q.ar sin i|i a v.r sin \|;f dm 

The integrals are evaluated over the volume of the rotor. In cylindrical 
coordinates, the differential mass is: 

dm * rdrdodz 

noting that : 


if ^ at ^ Q I J"[ ]dm « J" J J* [ Iprdxdrdc 


X r*0 cr*0 
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For brevity, let us perform the integration over disks (neglecting Che 
mass of Che rotor) only so that the dz just becomes "it," the axial thickness of 
the disk. So, now examine the integration with respect Co the face of the 
disk, e.g., dA ® rdrdo. For the case where the disk has no unbalance so that 
Che neutral axis, centerline, center of rotation, and centroidal axis are co- 
incident, it can be seen Chat:* 


g(x) 



2ir 

'l»\ /* /" 2 /cos(Sc+o)\ /0\ 

) dm * G(x) / I f ( Idrdor * I ) 

i|>/ J J \sin(nt+er)/ \o/ 

r 0=0 


similarly 


g(x) / r sin I* cos dm = G(x) If sin(Gc+o) coa(Gt+o) 

r 0 


drdo * 0 


and chat 


g(x) 


2ir 


f r^d« - 0U> f f P -3 

J \sin / J J \sin (Gt+o)/ 

r 0 


r drdo 


» G(x) 


pirr^ 


4 


vnien Che values of these integrals are substituted into the kinetic terms 
in the equations of motion, we obtain Che classical rotor equations. 

However, when Che disk center of gravity (CG) is no longer coincident 
with Che axis of .rotation, centerline, or rotor neutral axis, the integrals 
over Che face of the disk are no longer simple constants. For example: 

2 

hi y*r cos I* dm = hi(X(j)4x meG cos Gt =• Fyub* 


Funbt rotating unbalance, as "seen" by a stationary coordinate or 
support where me = unbalance. 


*Note Chat Che mode shapes are functions only of x, the axial distance. 


these functions be g(x) and 


/' 


g(x)dx 


G(x), 


% 


Let 
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When £he other sin cos ii, sin cos sin^tj cos^t integrals are 
evaluated over an unbalanced disk, so*^called secular terms are obtained 
terms multiplied by; 


sin at, cos Qt 

sin 2at, cos 2at 

If we recall that these integrals are coefficients of the generalized co“* 
ordinates, the equations of motion will have terms like these: 

q^sin at, q^cos at, q^cos2 at, q^sin2 at, q^sin2 at 

P.sin at, P.cos at, P.cos2 at, P.sin at, P,cos2 at 

r X ]L XX 

These terms are sometimes called parametric terms because they imply a 
harmonic variation in a parameter. In this case, these would be time^arying 
moments of inertia and mass moments* 

However, since these "unbalance” integrals are dependent only on the 
square of the unbalance, eccentricity, or noncongruency of the CG from the 
axis of rotation, the values are small and may be neglected. Of course, where 
the eccentric radius is no longer small, these parametric terms become impor- 
tant. 

Combining the terms from the casing, the connecting elements, and the 
rotor, the equations of motion of the entire turbine engine can be written 
down: 



X* 


(l)we assume that the rotor is symmetric so that the horizontal and vertical 
modes are the same and therefore all orthogonal. 
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( 1 ) 



22 


us: 



Recall that: M^j's are nonzero only for i * j; 
z*' - Zc(x*t) -^hicCx*) q£^ 

-yc(xtt) -2vi,(x*) Pi^ 

z* - Zr(xtc) - 23 

y* - 7r(xTt) - ^vi^(x*) ?i^ 


Therefore, the resulting matrix of equations will have an order equal 
to the number of all modes. ! 

The only excitations in this model (since we Ignored the parametric 
excitation) are those from the simple unbalances acting on the rotor: 



^m(x, 


o>e(xo)hi 


(Xf,)C0S 


Sit 



^^m(xo)e(xo)vi(xo)sin at 


Ah described in Vulum*! ILl, inlei I ci »?m •: loads .Much an rubs c«*in aJnn lin 
modeled and included in the TETRA computer code. The numerical solution for 
the transient response of a turbine engine with a sudden unbalance is also 
discussed in Volume III. Briefly, a central-difference integration approach 
is used to solve the governing equations. This consists of rewriting each 
ODE such that the generalized acceleration is on the left-hand side and the 
rest are to the right. The central-difference technique then calculates the 
current values of the acceleration, velocity, and displacement of each coord- 
inate in terms of three previous values . With these values, the value of the 
forces (right-hand side) is updated, and the calculation is repeated at each 
time step. 
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Examples of typical structural and connecting modules are given in 
Tables I through III- In TETRA, such modules and others are used to model 
an entire engine as well as its supporting Structure- 


Table I. Modal Module - 
Horizontal 



Static Structures. 
Plane 
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Table II. Modal Module Rotating Structures 


Horizontal Plane 



t 


^^11 

0 


0 


0 

K22 


0 


Mil 

0 

0 


0 

M22 

0 


Kll 

0 


0 


0 

^22 


0 



Horizontal 

Generalized Stiffness 



6 


Vertical 
Generalized Hass 



Vertical 

Generalized Stiffness 
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able II. Modal Module — Rotating Structure (Continued) 
Gyroscopic: Horizontal-Vertical Planes Coupling 




Table II. Modal Module - Rotating Structure (Concluded) • 
Unbalance Load - At Each Unbalance Location, x = x^. 




e h£ cos nt 
c h2 cos Qt 

c cos tXt 
e sin Ht 
e V2 sin Qt 

e Vq sin fit 
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Table III. Connecting Elements. 


Za = » £ hiaqi^ 


* * V 


'★ 

^h. q. 
1 j J-a J-a 


7 a = 7aCx,t) * 

* 3 ya , * 



-b » S ^ib <^ib 


7 b " Z 


- £ bb* qi^ 


9 b “ £ V-* Pi^ 


Linear Spring 




, q, - h. q. + h q. - h. q + "I 
la la lb lb ^a ^b *b J 


3 ^ ^ \ 

F = F rt ®F u. =u, k ju, p, -v Pi +U9P2 UoP2+*'> 

3o; V, 1, T, V 1/la lb -^b ^a ^a ” ^b b I 


Pja ®Pja 7a Ja Ja 

37 b * 

S F -r * F u. » -u. F 

Pjb yb 3pjb 7b Jb Jb 7a 
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Table III. Connecting Elements (Concluded). 
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3.0 DETAILED ANALYTICAL DEVELOPMENT OF THE TETRA PROGRAM 


3 . 1 SUMMARY 


This section contains the detailed application of the component element 
X3thcd ‘-.'7 Task II analyses for the transient dynamic response of a turbine 
engine. The component element approach results in a single independent vari- 
able, time, and reduces the number of differential equations required to obtain 
the essential dynamic characteristics of a turbine engine. The component ele- 
ments consist of elastic and rigid-body elements, described by generalized 
coordinates obtained by coordinate transformations, and physical connecting 
elements that model bearing/ frame springs and dampers, rotor-case rub springs, 
and gyroscopic crossaxis coupling effects. 

The generalized coordinates are based on the free-free modes and par- 
tially constrained modes associated with engine subsystem structures (rotor, 
casing, and pylon). The method extends conventional modal analysis procedure 
to acccount for physical damping, asymmetric stiffness terms, and rotor-case 
rubs including the effects of the dead band associated with the structural 
clearance. 

An efficient, numerical, time-integration scheme - the central finite-dif- 
ference method - is used to obtain the solution through an explicit, step-by- 
step calculation. 

A description of a computer code that was developed to implement the 
transient analysis method is also provided. 


3.2 INTRODUCTION 

This section is the detailed development of the application of the com- 
ponent element synthesis method described in the previous section for the 
transient dynamic analysis of turbine engine systems in three dimensions. In 
addition*, the structure and major program modules of a computer code that was 
developed to implement the transient analysis method is presented along with 
some computed results. 

The analysis method presented is conceptually similar to that described 
by Levy and Wilkinson (Reference 6) and clarified in Section 2.0 of this docu- 
ment. However, the method and the mathematical modeling approach described in 
this report have been specifically developed for the analysis of flexible tur- 
bine engine systems and include features needed to predict the response pro- 
duced by suddenly applied loads such as those caused by blade loss. 

The turbine engine is described by a reduced system of second-order dif- 
ferential equations, and a solution for the transient response is obtained 
through an explicit numerical-integration scheme; the central finite-differ- 
ence method. Global stiffness and damping matrices are not assembled, and only 
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the right“*hand side of the system of equations (i.e,, the forces) is updated 
at each time step. The differential equations are formulated in terms of 
generalized coordinates and modal elastic and rigid-body elements that repre- 
sent physical connections. In the former case, the elements are obtained 
through the coordinate transformations associated with the free-free modes 
and partially constrained modes computed for engine subsystem structures 
(rotors, case, and pylon). In the latter case, the elements describe the 
physical connections between the subsystem structures. These connections can 
be nonlinear and are defined with bearing/frame springs and dampers, engine 
support elements, link elements, rotor-case rub springs or stop elements, and 
gyroscopic crossaxis coupling elements. The rotor-case rub springs are used 
to model the additional load path, between the rotor 'and the case, that exists 
for the large rotor excursions caused by high rotor unbalance. For this 
ment, the effect of the force dead band associated with the structural clear- 
ance is included in the formulation for the effective restoring forces. 

The representation of the entire engine in terms of these components is 
achieved by a rational concatenation of the corresponding equations by appli- 
cation of Newton's Law or the principle of virtual work and the Lagrangian 
as described in Section 2.0. 


3.3 ANALYTICAL DEVELOPMENT . 

3.3.1 Governing Differential Equations 

The turbine engine system is described in terms of generalized coordi- 
nates and physical connecting elements. The generalized coordinates are ob- 
tained through coordinate transformations from the real space coordinates. 
Through the principle of virtual work, these generalized coordinates are the 
modal participation of specified engine subsystems. These modes and their 
generalized coordinates are associated with subsystems chat typically include 
the high and low pressure rotors, the engine case, and the pylon. Subsystem 
normal modes may be calculated from finite-element analyses or determined 
experimentally. Figure 3 shows a schematic of a subsystem and physical com- 
ponent* element definition for an engine system. The nonlinearities can be 
treated with the physical connecting elements, and the subsystem natural modes 
are used to define regions of the engine system which are expected to remain 
linear. As long as the nonlinearities are not severe enough to induce oper- 
ating regimes where the original modes are totally unrepresentative, the use 
of the natural modes is effective. 

Following the technique of the component element method, an ordinary dif- 
ferential equation is formulated for each generalized coordinate. These are 
written so that the acceleration is on the left-hand side and the rest are on 
the right, thus: 


“Kiqi = “Ci4i Fi(t) = Qi 


( 1 ) 
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The notation of Section 2,0 has been altered here to standard matrix no« 
tation in order to take advantage of matrix manipulation in the program coding. 
However, the principles are the same, 

qj are the generalized or modal coordinates that describe the subsystem 
modes. The q£ is the participation of the i-th mode of that subsystem* 
These are the principal dependent variables. 

m£ are the generalized or modal masses. 

k£ are the generalized stiffnesses that exist only for the distortion 
energy modes . 

* {*i}x {'liJ ^2) 

(K] is the physical stiffness matrix for the subsystem, and is the 

Eigen vector or mode shape of the i-th mode. 

are generalized damping coefficients proportionate to the modal 
masses 

C. 

1 Qpi mi 

where Qp. » the modal Q“f actor, and ui is the undamped natural frequency 
for the i-th mode. 

F£ are the generalized forces that are associated with applied physical 
forces. These physical forces include rotor unbalance forces, user- 
specified forces, crossaxis forces due to Coriolis acceleration, 
spring forces, damping forces, and rub forces, 

Qi are the total generalized forces which include the generalized re- 
storing forces - ^iqii the generalized forces - Ciq£, and the gen- 
eralized forces Fi associated with the applied physical forces. The 
total generalized forces are updated at each time step. 

It is useful to introduce the concept of a coupling ratio in formulating 
the generalized forces F£ in terms of the physical forces, 

^i ® ^r,d ®i,r,d 

The coupling ratio x d iorms the elements of {Q^} and is the deflec- 
tion at point r in the direction d for a unit value of the i-th generalized co- 
ordinate. The term F^ ^ defines the applied physical force at point r in the 
direction d. Then the*total generalized forces can be written as: 

Qi “ r d ^r,d ^i,r,d ^i^i ^i 4i (5) 
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and the equations of motion can be arranged into the form: 



1,2,3, 


( 6 ) 


3.3,2 Numerical Integration of the Governing Differential Equations 

The transient response solution for Equation 6 is obtained with an ex-“ 
plicit numerical’-integration scheme: the central finite-difference method. In 

this scheme, the accelerations of each of the generalized coordinates associ- 
ated with the component modes are computed at each time step using the current 
values of the applied forces. The accelerations are directly used to compute 
the future generalized displacements, and no iteration is required. The future 
displacements then allow the calculation of updated forces, and the procedure 
is repeated until the desired number of time-step solutions has been, obtained. 
As discussed in Reference 5, this concept has also been used in the prediction 
of crash response of vehicle structural systems. Global stiffness and damping 
matrices are never assembled. Local stiffness and damping matrices are de- 
fined, and it is only necessary to update those physical connecting element 
stiffness and damping matrices which are nonlinear. The forces are then de- 
termined by multiplication at the element level. This explicit scheme re- 
quires much less computer time and core storage than would be required in an 
implicit scheme where nonlinear system or global matrices must be redefined at 
each time step. 

Each of the equations of motion described by Equation 6 appears to be 
uncoupled. However, they are actually coupled through the physical constit- 
uent forces of F£ that relate to the displacement and velocity-dependent load 
paths between the modal subsystems that are used to define the linear regions 
in the engine system. 

In the central finite-difference method, the current total generalized 
forces in Equation 6 can be expressed- in terms of the current and previous 
values of the generalized displacements q£: 

Qi° - f (qi®. ql^ , ql2) (7) 

The superscripts 0, -1, -2 pertain to the current and two previous values 
of the generalized displacements (present value, value one time step earlier, 
and value two steps earlier). Then, the current generalized accelerations are 
obtained from: 

.. 0 Qi° (8) 

q . s . — 

1 m£ 

The following finite-difference equation is used to solve for the new or 
future generalized displacements: 

q.l = 2 qiO - qi'l + qiO (aT)2 (9) 
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The superscript 1 pertains to the future values of the generalized dis- 
placements • 


For convergent or stable solutions, the time step aT must be less than 


2 


where is the total effective stiffness that includes all product terms re- 
lating to qi in Equation 1. While not immediately apparent from inspection 
of Equation 1, the effects of q£ are included in F£ * It will be noted chat 
for nonlinear problems values are not known prior to the solution. In 
practice, aX should be made equal to about 1/40 of the smallest period of 
oscillation for the subsystem modes, assuming that there are no connection 
points, or: 



The initial conditions are accounted for by setting for the initial time 
interval : 

q 1 « q 0 + q 0 at 1 qO (aT) 2 
i i 1 7 

q_“l = q.O - q.O at + 1 q 0 (aT)2 (10) 

ITT 



^i 


.-2 . 2 qi-1 - q.O + qi^ ( T )2 


The latter condition arises from the assumption that the acceleration prior to 
time tQ is q£^« 

The major task in the implementation of Che method is the computation of 
the current total generalized forces From Equation 5, 


Qi° = ? d - ki qiO - Ci qiO 


( 11 ) 


The current generalized velocities are obtained from the following backward 
difference approximation: 


0 3 qi° - 4 q£~^ + "-“2 


^i 


( 12 ) 


2AT 
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Hence, in Equation (11), the force terms involving the generalized stiffness 
and damping can be readily defined • Conceptually, these force terms can be 
associated with scalar springs and dashpots connecting the generalized masses 
m£ directly to grounds The generalized force term 


includes coupling force terms between the generalized masses that are depen- 
dent on the current physical displacements and velocities of the physical con- 
necting elgments at the attachment points to the modal subsystems. In addi- 
tion, ®i,r,d includes nondisplacement and velocity-dependent terms that 

do not couple the generalized masses. These terms are associated with rotor 
unbalance and user— specif ied forces. 


The current physical displacements ^ and velocities ^ needed to 
establish the coupling forces are obtained from the summation of the modal 
participations of the generalized displacements. In other words, at each time 
step the modal displacements are transformed into real space (Equations 13 
and 14) so that the physical displacements can be obtained: 


2 n 

^r,d ® i ^i ®i,r,d 

2 

^r,d " i ^i^ ®i,r,d 


^r,d ^r,d physical displacements and velocities, 

at the point r in the direction d. 


(13) 

(14) 

respectively, 


The displacement and velocity-dependent terms of the forces associ- 

ated with a given connecting element are then determined at the local level by 
multiplying the connecting-element stiffness matrix into the physical dis- 
placements, the connecting-element damping matrix into the physical veloci- 
ties, and then adding the two results as depicted by Equation 15. It will be 
noted that in the case of the Coriolis acceleration forces only an effective 
local damping matrix exists for each station where this effect is included. 



2 2 

d r 


Xr,d Kr,d 





(15) 


Ky ^ and ^ are the current stiffness and damping coefficients for the 
connecting element. 


3.3.3 Modules and Structural Discretization 


The concept of the modal synthesis and the component element method 
naturally lead to modules or building blocks. These modules reflect physi- 
cal as well as mathematical components. With the first, the structure of 
the entire engine may- be constructed with physical subsystems as one may 
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sketch a schematic of the structure, and in the other a mathematical repre* 
sentation is built of the complete engine structure in terms of mathematical 
modules that mimic the response of the physical system. For example, it had 
been shown in Section 3,3,1, as well as in Section 2,0, that expansion of the 
physical displacement of each structural subsystem into its normal modes and 
application of Newton's Law or the Lagrangian results in modules of modal or 
generalized masses and stiffnesses analogous to the structural mass and stiff- 
ness of a simple spring-mass system. In addition, modal or generalized excit- 
ing forces can be constructed from the physical and external forces by multi- 
plying these with the modal displacement at the point. and line of action of 
the physical force. Rotation of a component also introduces gyroscopic mod- 
ules chat couple motions in two orthogonal planes. Similarly, the connecting 
forces have their own modules. In the following sections, the mathematical 
models of these modules are derived and written in forms amenable to computer 
coding. 


^•3.3.1 Modal Subsystems 

Modal representations obtained by transforming the degrees-of-freedom 
from real space to modal space are used to define the structural dynamic 
properties of the subsystems. This transformation can be made with the con- 
cept of the Lagrangian or virtual work principle as described in Section 2, 

The concept of using modal coordinates to reduce the dimensionality of the 
model and to clarify the dynamics of the model has historically been exploited 
in vibration analysis. Childs (Reference 7) and Lund (Reference 8) have re- 
ported success in using this approach for flexible rotor transient analysis. 


The equation describing the coordinate transformation for a given sub- 
system fellows: 


(x)- 


2 


<?i {<t^i) 


(16) 


{X} is the physical displacement vector whose elements are the ^ 
values 


qi are the generalized coordinates 


{^ 1 } is the mode shape vector whose elements are the values. 

The real-space formulations are based on beam-like, finite-element (dis- 
crete mass and stiffness) , subsystem models of the engine rotor(s) and case and 
a three-dimensional, finite-element model of the pylon, all described with re- 
spect to a ground-fixed coordinate system. Figure 4 shows typical subsystem 
normal modes and coupling ratios for a rotor finite-element model. The flex- 
ural vibration characteristics of equivalent nonrotating shaft(s) are used to 
define the rotor(s). The whirl phenomenon is addressed by using Euler's law 
for angular motions to establish crossaxis coupling forces that are dependent 
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on angular velocities in tvto planes and proportional to the polar inertia; 
these forces are applied as external forces, (contributors to the forces 
described in Equation 1) to the modal coordinates that represent the rotor 
-in two planes « 


The free“free undamped modes in both the vertical and the horizontal 
plane for the rotor(s) and the case are derived from planar finite-element 
models. In addition, a single-dimension (single degree-of freedom at each 
station) torsional model is also used to model the case. The rigid-body modes 
for the rotor (s) and case in the vertical and horizontal directions can be 
defined either with the **soft spring” rigid-body modes or with separate modal 
subsystems representing computed rigid-body modes based on mass properties 
and geometry. The torsional-direction, rigid-body modes for the case are 
represented either with the soft-spring, rigid-body modes or with a rigid- 
body modal subsystem. It will be noted that the case torsional direction is 
important for modeling the engine mounting system load paths. The fore-and- 
aft motions are represented with rigid-body modal subsystems for both the 
rotor(s) and the case. The pylon finite-element model is used to define a 
set of three-dimensional, cantilevered modes. Figure 5 shows the global 
coordinate system, direction numbers, and a typical engine mounting arrange- 
ment. Table IV identifies the physical degrees-of-freedom associated with 
each subsystem. Relative to defining the rigid-body mode shapes for separate 
rigid-body subsystems, it will be noted that, in general, an unconstrained 
rigid-body is free to move in space in six directions. Such rigid-body 
motion corresponds, in rectilinear coordinates, to translations in X, Y, and Z 


and rotations 
as follows: 


in 6 ys dg* equations describing these motions are 


My 


m 

0 

0 

0 

0 

0 


m 

0 

0 

0 

0 


m 

6 

0 

0 
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0 

h 

0 

0 


0 

0 

0 

u 


0 

0 

0 

0 
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Y 

•e 

z 

9x 

99 

6y 

99 


(17) 


Applied 

Physical 

Forces 


Inertia Matrix 


Physical Accelera- 
tion^ of the Center 
of Gravity 


These equations are included in the set defined by Equation 1. In this 
special case, corresponds to the physical mass properties m, I^, ly, or I 2 , 
q£ corresponds to the physical coordinates X, Y, Z, 8 y, or 62 of the 
center of gravity, and the k£ terms are equal to zero. Note that for a general 
(nonprincipal) set of axes, the inertia matrix of Equation 17 would include 
off-diagonal terms. For a three-dimensional body where the planes of symmetry 
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Thrust, Side, and Vertical 
at Forwa^rd Mount 


Side, Vertical, 
and Roll at 
Aft Mount 



Figure 5. T3rpical Mounting Arrangement and Coordinate System* 
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Table IV. Physical Global Degrees of Freedom and Direcrtion 
Numbers for the Subsystem* 



Global Degrees of 

Freedom 

Subsystems 

z 

9y 

B 


X 

0x 


n 

2 

3 

urn 

5 

6 

Vertical-Plane-Flexible 

Rotor(s) 

X 

X 





Horizontal-Plane-Flexible 

Rotor(s) 



X 

X 



Rigid-Body Rotor (s) 

X 

X 

X 

X 

X 


Vertical-Plane-Flexible 

Case 

X 

X 





Horizontal-Plane-Flexible 

Case 



X 

X 



Rigid-Body Case 

X 

X 

X 

X 

X 

X 

Torsional-Flexible Case 






X 

Three-Dimensional-Flexible 

Pylon 

X 


X 


X 























pass through the center of mass, there is no inertia coupling, and the inertia 
matrix corresponds to that of Equation 17, The Eigen vectors for the six 
rigid*“body modes are shown in Table V, 

In general the flexible, beam^like rotor and case finite-element models 
will be used to generate modal data for a single plane, say the vertical plane. 
These modal data can also be used to represent the horizontal plane. However, 
it will be noted that the sign of the modal data in either the Y or the Sg direc- 
tion must be changed in order to maintain a sign convention that is consis- 
tent with a right-hand, global-coordinate system. This is shown in Figure 6, 


3, 3,3,2 Applied Physical Forces 

In addition to the rotor unbalance and user-specified forces, applied 
physical forces include forces that are displacement and velocity dependent. 

The applied physical forces communicate with the generalized coordinates that 
represent the subsystems through. the coupling ratios, ^ The displace- 
ment and velocity-dependent forces act through the interconnections between 
the modal subsystems, and the load paths for these interconnections are defined 
with the physical connecting elements. 


3, 3, 3, 3 Physical Connecting Elements 

The force/displacement and force/velocity relationships for the physical 
connecting elements will be expressed in terms of Che stiffness and damping 
matrices shown in Equation 18, 

{F} =» - [K^] '{X} - [Cg] {X} (18) 

Tlie elements of {F} , {X}, and {X} define the terms F^^ X^ and X^ ^ 
respectively of Equation 15, > > » 

Equation 18 defines the physical forces exerted by the physical con- 
necting ^element on the modal subsystems to which it is connected. 


3, 3, 3, 4 Generalized Spring-Damper Element 

This element is associated with two physical points located at arbitrary 
locations in global space. Each of these points is assigned six degrees of 
freedom - three translational displacements (or velocities) and three rota- 
tional displacements (or velocities) , Thus, the dimension for either the 
stiffness or the damping matrix is equal to 12, 

This element can be used to represent simple, uncoupled, point springs 
and dampers that are useful in modeling rolling-element bearings. Among Che 
other load-path configurations this element can be used to model are engine 
frame/sump structures ,• hydrodynamic bearing definitions with direct and cross 
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However, Direction 0 is Opposite 

to Global Direction 0 » -0 

z 


Fieure 6. Vertical and Horizontal Planar Models. 



stiffness and damping, and crossaxis stiffness coupling arising from aerody«» 
naraic effects « A sketch of this element is shown in Figure 7, 

Typically, Che stiffness elements in [K^] are computed via a finite-ele- 
ment program, or with a closed-form solution for an idealized model, or are 
obtained from static or dynamic' testing. The damping elements in [C^] are 
computed if analytical expressions are available, or are obtained from dynamic 
testing, or are based on the assumption that the damping is proportional to 
Che stiffness. In this latter case [Ce]«[Kg], and the proportionality can be 
based on a specified percent of critical damping at a selected frequency as 
shown in Equation 19. 


[CJ 



(19) 




Specified Q-factor 


1 

2 C/Cc 


0 ) ® Selected frequency, radians/second. 


3*3. 3. 5 Space Link-Damper Element 

The SDiV-ze link-damper element sketched in Figure 8 is described by sixth- 
order stiffness and damping matrices. 
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-nm 


-Im 

nm 

ml 





m» — nl 

-Im 

-l2 

In 

ml 

l2 _ 



n, m, 

and 1 

are 

the direction 

cosines , and 

the column and 

row order correspond 


to directions 1, 3, and 5 at points I and J, respectively. 

A » cross-sectional area, in.^ 

E ^ Young *s modulus 
L ^ length 

The damping matrix [Cq] can be defined either in terms of translational 
(dashpot) damping directed along the axis of the link or by proportional 
damping. 





In the foraer case. 


C ^ specified scalar damping value, 

in. 

In the latter case, the form of the damping matrix is that shown in 
Equation 19, 


3.3. 3, 6 Rub Element 

The large rotor amplitudes associated with blade loss cause heavy rotor- 
case rubs that are usually accompanied by severe local damage to blades and 
case. It is to some degree fortunate that this local damage does occur, be- 
cause the loads acting on the bearings and frames are reduced by the action 
of the additional load path between the rotor and case chat is provided 
through Che mechanism of the heavy rubs. As a consequence of this action, 
the engine may be capable of withstanding, without catastrophic structural 
failure, the transient loading induced by the blade loss. 

The rub element allows the mathematical modeling of the nonlinear tip 
rub Chat includes the dead-hand displacement interval prior Co closure between 
the rotor and case. Upon closure, an equivalent linear spring, representing 
the local case distortion and blade compliance, is used to define the rotor- 
to-case load path. The net result is a bilinear spring with zero slope over 
the dead band and a finite slope over the region of interference. 


Figure 9 shows the relationships between the rub force and the position 
vectors of the rotor and case centers. Neglecting friction, the rub force 
can be defined as a force chat acts in the direction defined by the line of 
centers. The vector rub force acting on Che rotor can be written as; 


(1A| - c^) K - (V^ - V ) C for |A| > (22) 

r ! Aj o r c ' ' o 

® 0 for |A| < to 


Where K 


Radial spring rate representing the local case distortion (rep- 
resented by Che bulge shown in Figure 9 and the blade compliance. 


Zq ® Structural clearance 

^ Si (Dr «■ Dc) ® vector difference between Che position vectors of 
the rotor and case centers, 

7 * 77 “ ® unit vector in the direction of vector A, 
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At Time ■ t. Vectors, D and D 

c r 

Define Positions of Case and Rotor 
Centers. 

Vector Difference = a = D - D 

r c 

If IaI > e , A Rub Will Occur 
' ' o 

e = Structural Clearance 
0 

F =• Vector Rub » - — ^ (Ia| - e )K 

|Al ' ' ° 

K = Spring Rate Representing The Local Case Distortion and Blade 
Compliance 


Figure 9. Rub Force Model. 


and are the vector velocities of the rotor and case centers, re- 
spectively, 

C is the damping rate. 

In complex notation, A can be written as: 

A ® ** (yj. + j Z^) (y^ + j Z^) =* (y^- - y^) + j (Z^ - Z^-) 

Where j = unit vector in the z-direction. 

|A| = amplitude of A =* V^yr “ yc^^ *** ^ (23) 

* - [(yr - 7c> * j (2^ * Z^)] tl.O - eo/jAl ] K - (V^ - V^) C 

Fr * - (7r ” ^c^ (^-0 ^ So/Ul ] K - j (Z^ - Z^) U.O - zo/\^\ 1 ^ 

- C((y^ - y^)] - j C [Z^ - Z^)] 


Defining 

k ^ 1.0 


^o^ 

V^yr " 


K 


( 24 ) 


as an effective spring coefficient which is dependent on the initial struc- 
tural clearance and the rotor-case relative displacement and case-blade local 
spring rate. A stiffness matrix can be defined as follows: 



( 25 ) 


Where and F^y are the forces acting on the rotor center in the z and y 
directions , respectively. 


F ^2 p^y 3 J.Q che forces acting on the case center in the z and y di- 
rections , respectively. 

Zr, yr> Z^, y^ are the absolute displacements of the rotor and case 
center. 
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The stiffness matrix in Equation 25 is also shown in Figure 10 and repre-“ 
sents the [K^]* matrix in Equation 18. The damping forces can be written 
in matrix form as: 



(26) 


• « « • 

Zr, Yr, Zc, Yc are the absolute velocities of the rotor and case centers. The 
damping matrix in Equation 26 represents the [Cg]* matrix in Equation 18. 

It will be noted that the rub element defined by Equations 25 and 26 can be 
used to model rotor-to-rotor rubs as well as rotor-to~case rubs. In this case 
the inside rotor is identified with the subscript r, and the outside rotor is 
identified with subscript c. 


3.3.3. 7 Engine Support Element 

Real aircraft engine mounting systems are quite complex and must be mod- 
eled three-^dimensionally if accurate simulation is to be obtained. Figure 11 
shows an example of an aircraft engine mounting arrangement where three mount- 
ing planes are utilized. The forward mounting plane takes vertical, side, and 
axial loads. The midmounting plane takes side and roll loads. The aft mount- 
ing plane takes vertical loads only. 

The engine support element is a multipoint, multidirection element that 
provides the capability to model the complex load paths between the engine 
case and the pylon and thus allows for the direct modeling of actual engine 
mount structures. In addition, this element couples the flexible- and rigid- 
body-centerline modal subsystems that represent the engine case to the sup- 
port links that comprise the mounting system. Thus, this element is really a 
combined engine/support-case flexibility element that couples the beam-like 
casing model to the three-dimensional pylon subsystem (3-D). Figure 12 shows 
typical engine case-link attachment points. The loads at these attachment 
points induce frame/case distortions that increase the effective flexibility 
at the engine centerline. In this example, the side-direction flexibility at 
the forward mounting plane would be increased by the fan frame/case ovalization 
effect, and the vertical direction flexibility at the aft mounting plane would 
be increased by the turbine frame/case ovalization effect. 


Equations 25 and 26 represent the forces acting on the connecting sub- 
systems and already incorporate the minus signs shown in Equation 18. 
Hence, in this case, the [K^] and [C^] matrices shown in Equations 25 
and 26 should be multiplied by -1 to conform to the form of Equation 18. 
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Rub Element 



Rotor-Case Relative Displacement 
Vector Direction 


K - Spring Rate Representing 
Local Case Distortion 
and Blade Compliance 

£ = Structural Clearance 


Dead Band =* No Load Over 

Displacement Less 
Than £ 




Spring Rub Forces 

Acting on Rotor and 

Case When |a1 > e 
' ’ o 




£ 

O 


+ 





K 


Superscripts I and 3 Represent Directions Z and Y, 


Figure 10* Rub Element Model • 
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Vertical 


n Aircraft Engine 
ement • 


Vertical 



Turbine 

Mount 

Ring 


e local loads at the attachment points 
nting flexibility at the forward plane, 
nd at the aft plane, in the vertical. 

e/Case Ovalization Effects Increase the 
exibilitv at the Mounting Planes. 






Figure 13 shows the attachment and centerline points in global space 
and Che load directions for Che two mounting planes corresponding to Che typi« 
cal mounting arrangement shown in Figure 5, At the forward mounting plane, 
points L, M, and N are. the attachment points on the pylon, and point I is at 
Che engine centerline. The two attachment points on the engine case are iden- 
tified by J and K. In this case, the engine support element is represented 
with the stiffness matrix relating Co Che points I, L, M, and N. The dis- 
placement terms associated with points J and K are reduced out. Point I has 
six degrees of freedom, and points L, M, and N each have three degrees of 
freedom. It is necessary to eliminate points J and K because the physical 
connecting elements are only used to define the load paths between modal sub- 
systems or between modal subsystems and ground. 

The model definition for the portion of Che engine support element Chat 
pertains to the case flexibility (flexibility between the engine centerline 
and Che case-attachment points) is shown in Figure 14. In this representation 
all of the case flexibility is lumped at the case-attachment points (Points J 
and K) , and rigid members are used to transfer the loads to the centerline 
point (Point I). The case flexibility is described by three spring rates K^, 
K^, and K^. These rates are reciprocals ’of the case flexibility values 
obtained by either calculation or measurement. In the former case, a finite- 
element model of the frame is typically used to define the flexibility values. 
For example, Ky would be defined by applying an n ® 1 shear flow and restrain- 
ing the frame model at Che case-attachment point locations. The vertical dis- 
placement of a horizontal diameter would then be used to establish Che spring 
rate Ky. Figure 15 shows the twelfth-order stiffness matrix that describes 
the case flexibility. The terms a, b, c, d, e, and f are multipliers used to 
proportion the spring rates, and Che associations between the spring races Ky, 
K^, and the lumped spring location points J and K are shown below. 

Case Distortion Spring Rates 

Ky =• Vertical Spring Rate 
» Horizontal Spring Rate 

K^ =* Axial Spring Rate 


Multipliers 


Spring Races 
(Ib/in.) 

Point 

J 

K 

Kv 

a 

C 

Kh 

b 

d 

Ka 

e 

f 


In conjunction with this scheme, the following restraints must be 
followed. 

a ^ c * i.O 


b d « 1.0 
e*^f— 1.0 

For symmetrical mounting, a®b-c®d = e=®f®0.5. 




Combined Engine Support-Link Element 
Multidirection, Multipoint Element 
Represents Case Distortion and Links 
Flexibility . 
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Distortion and Links Flexibility. 




Forward 


Direction 

5 


Engine 
Frame/Case 


(Forward, Looking Aft) 


Numbers 

This Spring Element 
Couples, the Flexible 
and Rigid-Body 
Centerline Subsys- 
tems That Represent 
the Engine Case to 
the Case- Attachment 
Points 

The Locations of the 
Case-Attachment 
Points Are Arbitrary 
At the Mounting 
Station 


Case-Attachment Pointy 

"Rigid Frame** Transfer^^ 
the Centerline Loads to 
the Lumped Springs at 
the Case-Attachment 
Points at J and K 



Case- Attachment Point 


© 

Engine Centerline 


^ Forward 


The Frame /Case Flexibility is Lumped at the 
Springs Shown at Points J and K 


The Lengths 
and d Define*^ the 
Case-At t achment 
Points Relative to the 
Engine Centerline 
Located at Point I 


Figure 14. Case Flexibility Represented by Multipoint/Direction Physical 
Spring Element (Which Defines a portion of the Engine Support 
Element) . 
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The engine support element is a variable-geometry element formed by com- 
bining Che spring element stiffness matrix [Kf] of Figure 15 with Che stiff- 
ness matrices of the connecting link elements. The stiffness matrices for the 
individual link elements have Che same form as Che stiffness matrix shown in 
Equation 20. 


[KJ 

[Kel 

[Kf] 

[Ki] 

[KcJ 

tKeJ 


[Kf] + [Kt^] 


= r^i_i.5i2 


[Kii - Ki 2 K22"i K21] 

Spring element stiffness matrix representing the case 
Stiffness matrix for the links 

Total unreduced stiffness matrix for the assembly 
Reduced stiffness matrix for the engine support element 


(27) 

(28) 


The matrix [K^,] in Equation 27 represents the total unreduced stiffness 
matrix and includes the case Points J and K* Equation 28 represents the 
matrix reduction that is performed to eliminate Points J and K. Figure 13 
shows two examples of engine support element configurations. It will be 
noted that there are no local moment load paths at Points J and K or at the 
link attachments at the pylon. 


The damping matrix [C^] for the engine support element is defined in 
terms of proportional damping and has the form of the damping matrix shown 
in Equation 19. 


3. 3.3.8 Gyro Element 


The gyro element models the crossaxis forces due to Coriolis accelera- 
tion and addresses general whirl motion defined by the response in two planes. 
These cfossaxis forces are mathematically treated as "right-hand-side forces," 
or externally applied forces, and couple the vertical and horizontal-plane 
rotor subsystem models. Figure 16 shows the damping matrix used to define the 
gyro forces acting on the connecting subsystems. This matrix represents -1 
times the [Cg] matrix in Equation 18 and is derived from Euler’s equations for 
rotational motion. For a spinning disk on a whirling shaft, Euler’s equations 
of motion in a fixed frame can be written as: 



( 29 ) 
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Figure 15, Combined In-Plane and Out-of-Plane Stiffness Matrix for the 
Physical Spring Element Shown in Figure 14. 
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Figure 16 • Gyro Element. 


My and M 2 are Che applied momenCs , 



The inertia matrix that defines Che uncoupled rotary 
inertia moments acting in the two planes. 



The damping matrix that defines the coupled-velocity 
dependent gyroscopic moments. 


The effects of Che uncoupled rotary inertia are implicitly included in 
the rotor subsystem modal data, and the velocity-dependent moments are treated 
as the applied physical forces shown in Figure 16. 


3.3 .3.9 Unbalance Forces 


At a given rotor mass station, the unbalance forces are related to the 
radial separation of Che geometric and mass centers associated with Che speci- 
fied unbalance. The unbalance forces are inertia-type forces of a transla- 
tional nature contrasted to Che rotary and gyroscopic forces which are of a 
rotational nature. This allows us to make use of an important theorem in ana- 
lytical mechanics which states that translational motions of the center of 
gravity of a rigid body are not coupled by inertia forces to Che rotations 
of the body about its center of gravity. This theorem is attributed to Euler. 
Consequently, we may separate the translational and rotational inertia forces 
into two groups and treat them independently. The effects of Che rotational 
inertia forces are incorporated via the subsystem modal data and through Che 
gyro elements as mentioned previously. The unbalance forces associated with 
Che translational motions of the center of gravity are identical in form to 
those acting on a mass particle. 

Figure 17 will be used in the development of the acceleration forces due 
to disk translation. 

Defining: c » Radial separation of geometric center and the mass center of 

^ Che disk. 

^2 ® Unit vector in Z direction of fixed reference frame. 

Hy ® Unit vector in Y direction of fixed reference frame. 

The vector velocity of the mass center can be written: 

Velocity of mass center « ® ™ (Y + e cos wt) n ~ (Z + e sin wt) n 

T dc “ y ^ z 

Replace ut with 9 V_ = (Y + e cos e) n + (2 + c sin 9) n 

T dc “ ^ dt “ ^ 
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" Su 9 9^ riy + (Z + Su 9 cos 0) rij 

» vector acceleration of mass center. 


d * " d * * 

At » (Y - 0 sin 9) — (z + 0 cos 0) 

At = (Y - Cy 0 sin 0 - ey 0^ cos 0) riy 
CZ ^ Ey 0 cos 0 * s ti 9^ sm 0) Tit 

• • V« 

The Y and Z acceleration terms are already included in the modal subsys- 
tem representation of the rotor. The remaining accelerations are as follows: 


Acceleration Term 


Z Direction Y Direction 


Centripetal Acceleration 9^ sin 9 9^^ cos 9 


Angular Acceleration Su 9 cos 9 -Cu 9 sin 9 

The fixed frame reversed effective or D* Alembert forces applied to the 
mass stations of the modal subsystems representing unbalanced rotors are of 
the following form: 


Z Direction (Vertical-Plane Subsystem) 

Fg .« M 9-^ sin 9 - M e^ 9 cos 9 = Uui^ sin wt - U a cos a>C 

Y Direction (Horizontal-Plane Subsystem) 

Fy ® M 9^ cos 9 + M 9 sin 0 « U cos ut + U a sin ut 

M ® disk mass 
U » M €^ ® mass unbalance 

3 e -T” (cii) ® angular acceleration 
^t 

Including the effect of unbalance phase, as shown in Figure 18, the 
applied physical forces in the fixed frame due to unbalance at Point r are 
written as follows : 
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Fr ^2 = force at point r in Z direction 

“ Uy 0 )^ sin (wt + -*• a cos (cat <j)y) (30) 

Fr,y ** force at point r in Y direction 

® cos (oit a sin (ojt + (31) 


= Unbalance Vector- 



Figure 18. Effect of Applied Physical Forces: in the Fixed Frame, 


In order to define the unbalance forces, the current angular velocity 9 
and angular acceleration of the unbalanced rotor, must be known. In addi- 
tion, the current angular velocity of the rotor is> needed to calculate the 
gyro moments. The angular, position, velocity, and. acceleration of the unbal- 
anced rotor are treated independent of the time-transient, systemrresponse 
values. The speed history of the rotor is defined in terms of an accel/decel 
rate and an ending time for a series of speed segments. 

0Q » currrent angular position of rotor. 

9i =« future angular position at time. t]^. 

tg. ® current time 

At - time step 

ti » future time - tQ ^ 
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® angular acceleration for time interval At 
* future angular velocity 

• 1 2 

“ ®0 ®0 (ac) (32) 

®l “ ®0 (33) 


3 •3.3,10 Subsystem Internal Loads 

The internal loads for a given subsystem are determined by superposing 
the effects of the applied physical forces that are displacement or velocity 
dependent and the physical forces derived from the modal contributions* In 
the former case, the physical forces are related to changes in the boundary 
conditions to the original modal subsystems; in the latter case, Che physical 
forces are associated with transformations from modal to real space. 

Total internal loads: 

^r,d ^ ^r,d^ 

where ^r,d “ contribution of applied physical force 
=* contribution of modal forces 

^r,d =» 2 ^^i,r,d^ 

i 

where ^ represents the force at Point r in the Direction d for a unit 
value of*tlie i-th generalized coordinate. 

It will be noted that the values of the generalized coordinates are 
dependent on the applied physical forces that are displacement and velocity 
dependent as well as the applied physical forces that are not dependent on 
displacement and velocity. 

Examples of the internal bearing forces are given in Figures 41 and 42 in 
Section 3.5. 


3.4 STRUCTURE OF THE COMPUTER PROGRAM 


The TETRA computer program automatically .combines subsystem modal data 
and physical connecting element data to form the governing differential equa*** 
tions described in Equation 6 for each of the generalized coordinates and per- 
forms the numerical integration. The major subroutines where the transforma- 
tions and the numerical integration of the generalized differential equations 
are performed are identified in Figure 19. 
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CITRRT 

Performs Transformation 
From Generalized To 
Real Space Needed 
Co Establish 
Physical Forces 





0 

The Cask of establishing Che current tocai generalized forces Qf acting 
on each of the generalized coordinates is accomplished in the subroutine 
GEN. The current physical displacements and velocities, needed to establish 
Che physical coupling forces used by GEN in the formation of the generalized 
forces, are determined in subroutine CURRT which transforms Che modal dis- 
placements qi into real space at each time step. It will be noted chat, in 
computing the generalized forces from the physical displacements and veloci- 
ties, a transformation from real space to generalized space is performed by 
subroutine GEN for each time step. The generalized forces also include the 
effects of applied physical forces which are not displacement or velocity 
dependent, such as unbalance forces. The generalized forces and the genera- 
lized masses are utilized in subroutine TILOOP to compute the generalized 
accelerations using Equation 18 and to compute the future values of the gen- 
eralized displacements using the central finite-difference method as shown 
in Equation 9. 

The modal subsystem numbers and the types and the regions they represent 
are identified in Table VI. The number of directions (physical degrees of 
freedom, DOF) and the direction numbers associated with these degrees of free- 
dom are listed in Table VII. 


Table VI. Modal Subsystems. 


No. 

Subsystem 



Type 

1 

Vertical Plane Rotor - 

1 

Model 

Flexible-Planar 

2 

Horizontal Plane Rotor 

- 

1 Model 

Flexib le-P lanar 

3 ■ 

3-D Rigid Body Rotor - 

1 

Model 

Rigid with 5 DOF 

,4 

Vertical Plane Rotor - 

2 

Model 

Flexib le-P lanar 

' 5 

Horizontal Plane Rotor 

- 

2 Model 

FI ex ib 1 e-P 1 anar 

6 

3-D Rigid Body Rotor - 

2 

Model 

Rigid with 5 DOF 

7 

Vertical Plane Case Model 

Flexib le-P lanar 

3 

Horizontal Plane Case Model 

Flexible-Planar 

9 

3-D Rigid Body Case Model 

Rigid with. 6 DOF 

10 

Torsional Case Model 



Flexible-Twist 

11 

3-D Flexible Pylon (Wing) 

Model 

Flexib le-3-D 
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Table VII, Subsystem Identification 


Subsystem Model 


Subsystem 



Description 


Rotor-1 Vertical Plane 
Flexible Model 

Rotor-1 Horizontal Plane 
Flexible Model 

Elotor-1 Rigid Body 
Model 


Rotor-2 Vertical Plane 
Flexible Model 

Rotor-2 Horizontal Plane 
Flexible Model 

Rotor-2 Rigid Body 
Model 


Case Vertical Plane 
Flexible Model 

Case Horizontal Plane 
Flexible Model 

Case Rigid Body Model 


Case-Torsional Flexible Model 
Pylon 3-D Flexible Model 


Directions 

Number 

Direction Numbers 

2 

K»1 


K=2 

2 

K-3 


K»4 

5 

K»1 


K=2 1 


K»3 




K»5 

2 

K-1 


K*2 

2 

K=>3 


K*4 

5 

K»1 


K=2 


K»3 


K«4 


K=5 ! 

2 

K»1 


K-2 

2 

K=3 


K=4 

6 

K=1 K=*4 


K»2 K»5 


K=3 K=6 

1 

K=6 

3 

K=l K=5 ^ 


K«3 












Eleven triple-subscripted arrays are used to store the subsystem mode 
shape data: 

Si (I,J,K) » Mode shape modal displacements for the 
i-th modal subsystem. 

i * 1, 2, 3, . . . .,11 (Table VI) 

I identifies the generalized coordinate 

J identifies the location or point on the mode shape 

K identifies the direction (Table VII) 

The Si array data are used to define the coupling ratios Bi ,^,4 needed to 
perform the transformations between real and generalized space. 

The physical displacement and velocity data for the current time step are 
stored in the arrays X(I,J) and VEL(I,J), respectively: 

I * point number 

J » direction ntanber (1-6) 

The generalized coordinate values are stored in the array Z(I,J): 

Z(I,J) » Present and previous values for the generalized coordinates 
I » Generalized coordinate number 

J ® 1 for current time (0) , 2 for one previous time step 

(-1), 3 for two previous time steps (-2). 

The generalized masses, stiffnesses, and damping coefficients for the 
governing differential equations are stored in the arrays ZM, ZK, and ZC: 

ZM(I) » . Generalized masses 

ZK(I) ® Generalized stiffnesses 

ZC(I) ® Generalized damping coefficients 

The physical coupling forces used in the formation of the generalized 

forces are computed by multiplying the physical displacements and velocities 
into the connecting element stiffness and damping matrices, respectively, and 
then adding the results. The current values for the physical connecting 
element data are stored in Che arrays, AKE^ and ADE^. 

AJ^i(I,J,K) ^ Stiffness array for the i-Ch physical connecting 
element type. 

ADEi(l,J,K) » Damping array for the i-th physical connecting 
element type . 
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i ® Is = Identifies the physical connecting element. type« For 

example, the generalized spring-damper element, the 
rub element, etc, 

I Identifies the physical connecting element number, 

J » Identifies the row number associated- with a force 
for a given point and direction, 

K - Identifies the column number associated with a displace- 
ment for a given point and direction. 

The information needed to identify the locations of the AKE and ADE 
array elements is provided in the following arrays: 

ICOMPJ (I) ® Number of points for Element I, 

ICOMPE (I,J) ® The global point numbers for Element I associated with 
the J-th point, 

ICOMPD (I,J) » The number of directions for Element I* associated with 

the J-th point, 

ICOMPNd, J,K) » The direction numbers for Element T associated with the 
J-th point and the K-th direction, 

3,5 TURBINE ENGINE SIMULATION 

The results obtained with a demonstrator model based on a VAST simulation 
of the single-spool engine system in Figure 20 will now. be presented. The 
straight lines shown for the rotor and case in this figure represent segment 
beam-centerline elements that include both stiffness and mass- properties. 
Free-free modal data for the rotor and case subsystems along with physical, 
spring-damper connecting element data were input to the TETRA time-transient 
computer program, and solutions were obtained for various rotor speed and 
rotor unbalance conditions. 

Modal data based on planar finite-element models of the total system 
shown in Figure 20 and on the free-free rotor and casing subsystems were com- 
puted with the VAST program. In the former case,, a frequency-domain modal 
analysis was used to establish the steady-state frequency response of the 
total system: In the latter case, the TETRA program was used to synthesize the 

free-free modal data and the physical connecting element data and predict the 
time-domain response for the combined system. 

Spectral analysis and a comparison of the TETRA and frequency-domain 
solutions indicated that TETRA had correctly synthesized the modal and physi- 
cal data to predict the time-transient response for the combined system. 
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Fan Turbine 

(Gyro Station) (Gyro Station) 


Physical Connecting Elements 

Spring 1, 0.5 x 10® Ib/in. 
Spring 2, 0.5 x 10® Ib/in. 
Spring 3, 1.0 x 10® Ib/in. 
Spring 4, 1.0 x 10® Ib/in. 


Point 1 is the Global Origin 

(X = Y * Z a 0) for the TETRA Model 

Physical Point Numbers 1, 2, 37 
etc. Identified 


Figure 20. Finite-Element Model for Demonstrator, 40 DOF. 
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The frequencies for Che free-free subsystem modes used in Che TETRA 
analysis are shown in Table VIII, Eight generalized coordinates, correspond- 
ing to these frequencies, were used to define the time response for the total 
system which has 40 degrees of freedom. Therefore, the dimensionality of Che 
model has been reduced by 80%. 


Table VIII. Frequencies^ for the Free-free Modal 

Subsystems Used in the TETRA Analysis. 


Rotor 


Case 


cpm 

Hz 

cpm 

Hz 

99.3 

1.655 

125.1 

2.085 

120.8 

2,013 

159 

2.65 

10065 

167.75 

30367 

510.62 

16851 

280.85 



26425 

440.42 




"^Free-free modal data obtained with the following boundary 
conditions. The rotor was supported at each end with 300 
lb/ in. spring, and the case was supported at each end with 
a 100 lb/ in. spring. 


Inspection of Table VIII shows that the smallest period of oscillation 
for the subsystem modes is equal to 1.0/510.62 » 0.001958 second. Account- 
ing for the stiffening effect of Che physical connecting springs, a time step 
value tKat is considerably less than 0.001958 second must be used for the 
numerical integration Chat will be performed with the TETRA program. In 
practice, AT should be made equal to about 1/40 of the smallest period of 
oscillation for the subsystem modes: 

i 3t 0.001958 =» 4.895 x 10~5 

40 

On this basis, a time step of aT ® 50 microseconds was selected. 

The free-free modal data were used in the TETRA program to represent both 
the vertical and the horizontal plane, and uncoupled springs (in the Z and Y 
directions) were used to connect the modal subsystems to each other and to 
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the ground to model the configuration shown in Figure 20. The rigid-body 
modes for the rotor and the case were approximated with the soft-spring, 
rigid-body modes at 1.655, 2.013, 2.085, and 2,65 Hz. In addition, they were 
defined in subsequent analyses with the true representations based on the mass 
properties. The results obtained with the two approaches were in good agree- 
ment. A modal Q-factor equal to 15 (3.33% critical damping) was used to rep- 
resent the damping of each of the casing subsystem modes. The rotor subsystem 
.nodes were undamped. Proportional damping based on a physical Q-factor equal 
to 15, a selected frequency of 3624 rpm, and corresponding to one of the crit- 
ical speeds was used to represent the damping for the connecting spring 
elements. 

For constant-speed running, the following time-step and total-time param- 
eters were used in the TETRA time-integration analysis: 

AT « 50 lAsec 

Number of time integration points « 10,240. 

Time frame * 10,240 x 50 x 10~^ =* 0.512 seconds. 

Figure 21 shows the TETKA computed displacement-time history at the 
middle of the rotor (Point 38) in the vertical direction corresponding to a 
sudden 100 g-in. fan unbalance at a constant 12,000 rpm speed. The number 
of computed values was decimated so that 1024 points are shown, and the time 
increment between points is equal to 5 microseconds. The 12,000-rpm speed 
corresponds to supercritical-speed operation relative to- a rotor-dominated 
mode computed at 9908 rpm for the total system. Overshoot is clearly in evi- 
dence for the transient response shown in Figure 21. The fast Fourier trans- 
form of this time response, shown in Figure 22, indicates that four modes are 
contributing to the response. The peak at 200 Hz (12,000 rpm) is associated 
with the driving force. The VAST predicted mode shapes for the total system 
sketched in Figure 23 show that the modes at 3,617, 9,908, 13,983, and 26,473 
rpm should be responsive at the middle of the rotor for fan unbalance. This 
indicates that TETRA has correctly synthesized the modal and physical data to 
predict the time- transient response for the total system. Figure 24 shows the 
displacement-time history at the middle of the rotor in the vertical direction 
corresponding to a sudden 100 g-in. fan unbalance at 9908-rpm critical speed 
operation. In this case there is no overshoot, and the transient response 
builds up to the steady-state level. This response is characteristic of oper- 
ation at a critical speed. 

The^ results presented in Figures 21 and 24 do not include gyroscopic 
stiffening. Figure 25 shows the off-resonant response in two planes at the 
fan that reflects the effects of gyroscopic stiffening at both the fan and 
turbine. Figure 26 shows orbit plots of the data presented in Figure 25. 
Inspection of Figure 26 shows that the initial response is a noncircular 
whirl. This noncircular whirl transitions into a circular whirl after steady- 
state conditions are reached. Figures 27 through 32 show response plots at 
the fan and case that include both the effects of gyroscopic stiffening and 
the rub load path. Figures 33 through 36 show the steady-state frequency 
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Figure 21. Demonstrator Model 100 g-ln. Sudden Unbalance at the Fan at 12,000 rpm. 




Driving Force 


Response 

Point Rotor Speed 



Figure 23. Total System - Rotor Mode Shapes 
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Displacement in the Z Direction 
at the Middle of the Rotor , in. 


A T = 5 Msec, 1024 Points 



Figure 24 • Demonstrator Model 100 g-in* Sudden Unbalance at the Fan at 
9908 rpitte 
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(T =* 0.105 to 
T = 0.3 sec) 


Figure 26. Orbit Plots at the Fan for NASA Demonstrator Model with Gyro. 
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Figure 27. Response In the Vertical Direction at the Case and 
Fan Rotor for the NASA Demonstrator Model for 5000 
g-in. Sudden Fan Unbalance at 3000 rpm (No Rub). 
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Vertical Displacement, in. 


Shown to Time 0.05 seconds^ Radial-Displacement Dead Band 
Exceeds the Fan Rotor-Case Relative Displacement, 



-0.08 “0.06 “0.04 “0.02 0 0.02 0.04 0.06 0.08 

Horizontal Displacement , in. 

Figure 29, Rotor and Case Orbit Plot for the NASA Demonstrator Model 
for 5000 g“in. Sudden Fan Unbalance at 3000 rpm (No Rub) . 
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Vertical Displacement , in. 


Shown to Time 0.05 seconds; iO-mil , Radial-Displacement Dead 
Band and 1 x 10® Ib/in. Rub Spring at the Fan Rotor Case. 



-0.08 -0.06 -0.04 -0.02 0 0.02 0.04 0.06 0.08 

Horizontal Displacement, in. 

Figure 30. Rotor and Case Orbit Plot for the NASA Demonstrator 
Model for SOOO g-in. Sudden Fan Unbalance at 3000 
rpm (With Rub) . 
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Fan Displa< 



Figure 33. NASA Demonstrator Steady-State Frequency Response for 
100 g-in. Fan Unbalance (No Gyro) . 




Middle of Rotor Displacement/Peak Middle of Rotor Displacement 



0 3000 6000 9000 12,000 15,000 

Speed, rpm 

Figure 34* NASA Demonstrator Steady-State Frequency Response for 
100 g-in* Pan Unbalance (No Gyro) . 
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Fan Displacement /Peak Fan Displacement 



Figure 35. NASA Demonstrator Steady-State Frequency Response for 
100 g-in. Fan Unbalance (Gyro at Fan and Turbine) . 
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Figure 36. NASA Demonstrator Steady-State Frequency Response for 
100 g-in. Fan Unbalance (Gyro at Fan and Turbine). 


response for the total system with and without gyroscopic stiffening at the 
fnn and turbine locatinuM. (»f * Im* r iiiM* n-Mpon/M*, nucf* m t y-nl ati* 

conditions have been reached with the frequency response, shows good agreement* 

Figures 37 through 42 show the results of analyses for a 2000 rpm/sec 
accel rate that reflects use of the TETRA restart option* The restart option 
allows the utilization of the results obtained from a previous analysis to 
continue without the introduction of pseudotransients . Each time output is 
printed, the TETRA program writes in a file the values for each of the gener«» 
alized coordinates for the last time step, for one time step earlier, and for 
two time steps earlier* This information is utilized to provide the initial 
data needed to continue the analysis* The process can be repeated any number 
of times. 


87 




Time, seconds 

Figure 37. Response in the Vertical Direction of the Case and 
Pan Rotor for the NASA Demonstrator Model for 1000 
to 3000 rpm Accel Segment; Radial -Displacement Dead 
Band Exceeds the Pan-Case Relative Displacement 
(No Rub). 
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Figure 39. Response in the Vertical Direction at the Case and 
Fan Rotor tor the NASA Demonstrator Model for 1000 
to 3000 rpm Accel Segment; 10-mil Radial-Displace- 
ment Dead Band and 1 x 10® Ib/in. Rub Spring at the 
Fan Rotor-Case (With Rub). 
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Fan Speed, rpm Force, Ibf 
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Figure 41. Force Exerted by the Forward Frame /Bearing (Spring 3) 
on the Rotor at Point 5 in the Vertical Direction for 
the NASA Demonstrator Model for 3000 to 5000 rpm Accel 
Segment; Radial-Displacement Dead Band Exceeds the Fan 
Rotor-Case Relative Displacement (No Rub). 



Fan Speed, rpra Force, Ibf 


5100 g-in. Fan Unbalance 


4040 rpm 



Time, seconds 


Figure 42. Force Exerted by the Forward Frame/Bearing (Spring 3) 
on tbe Rotor at Point 5 in the Vertical Direction for 
the NASA Demonstrator Model for 3000 to 5000 rpm Accel 
Segment; 10-mil Radial-Displacement Dead Band and 1 x 
10® Ib/in. Rub Spring at the Fan Rotor Case (With Rub). 




The computing strategy used to generate the response data shown in 
Figures 37 through 42 was as follows* The maximum speed for the accel was 
selected as 5000 rpm to encompass the two gyro-stiffened critical speeds 
located at 3292 rpm and 3624 rpm for steady-state operation* Modal Q-factor 
values of 15 were used for each of the casing subsystem modes, and the rotor 
subsystem modes were undamped* Physical damping for the connecting springs 
based on a Q-factor of 15 and a 60.4-Hz frequency, corresponding to the 3624 
rpm mode, was used* 

The analysis was accomplished with three restart segments for the follow- 
ing cases: 

• Large fan-case radial clearance (200 mils) 

• Small fan-case radial clearance (10 mils) 


Initial constant-speed runs were made for 0,520 second of simulated time with 
the rotor speed set at 1000 rpm and with 100 g-in, fan unbalance to establish 
steady-state conditions. For the 2000 rpm/ sec accel rate, the time required 
to accel from 1000 to 5000 rpm is 


(5000 - 1000) rpm 
2000 

sec 


2 seconds. 


With a integration time step of aT = 50 ^sec, the number of time steps required 
is equal to 


50 


2 seconds 
microseconds 
time step 


40,000 time steps 


For both cases, two restart segments of 20,000 integration points each were 
used to compute the response for the 2-second accel from 1000 to 5000 rpm. 
For plotting purposes, the computed data were decimated so that each of the 
restart segments plots represents 2000 points with a time increment between 
points equal to 5 microseconds. This means that there are at least: 


— X ^ ® 24 

5000 5 X lO"*^ 

plotted points per cycle at the highest rotor speed (5000 rpm). Although the 
plots are labeled 0 to 1.0 second for the accel from 1000 to 3000 rpm, the 0 
time actually corresponds to 0.52 second. At a time 0.6 second (labeled as 
0.6 - 0.52 » 0,08 second), a large unbalance increment of 5000 g-in. is added 

to represent sudden unbalance. It will be noted that the apparent critical 
speeds for the large-clearance case (i.e., no rub) occur at higher rotor 
speeds than the steady-state values. The 3530 rpm speed corresponds to the 
3292 rpm speed, and 3960 rpm corresponds to the 3624 rpm speed. The shifting 
of the peak-response speeds is caused by the accel rata, and this phenomenon 
is discussed in Reference 9. The fan rub present for the small-clearance case 
causes a stiffening effect that shifts the peak response to a higher speed. 
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4.0 APPLICATION OF TETRA TO BLADE LOSS TRANSIENT RESPONSE ANALYSIS 
OF TWO TEST VEHICLES AND EXPERIMENTAL VERIFICATION 


4.1 SUMMARY 


This section discusses the results of Che application of TETRA to two 
blade-out test vehicles and the test data obtained from these tests. Included 
are the descriptions of these two vehicles, the comparison of dynamic responses 
(predicted versus test), and an assessment of the TETRA computer program. 


4.2 INTRODUCTION 


As Galileo and Sir Francis Bacon said centuries ago, a theory stands or 
falls by its ability to predict experimental results. The test of usefulness 
or validity of any calculation method, therefore, must be twofold: (1) pre- 

diction of the result of an experiment and (2) experimental verification by 
comparing predicted with measured values. 

Because the TETRA program had been checked during its development with 
simple, degenerate cases wherein known closed-form solutions are readily avail- 
able, this section will concern itself with the response of complex structures. 
It is part of Task II to verify the computer program developed in this task 
by applying it to the calculation of the transient response of a test vehicle 
and an engine where dynamic test data had been obtained during blade-out tests. 

This section describes two experimental vehicles and presents the dynamic 
response calculated by TETRA under the same transient, blade-loss loading of 
the tests. Loads and displacements are calculated at points corresponding to 
the locations of instrumentation in the test vehicles. 

Task III effort comprises obtaining and reducing the test data from these 
two vehicles and preparing them in a form to permit comparison with the pre- 
dicted response. 

Finally, the effectiveness of the TETRA program is assessed in Section 
5.0. Refinements to improve its accuracy and to reduce its shortcomings, as 
well as extensions to improve its versatility and increase its scope, will be 
recommended. 

Results from two separate blade-out tests are used to provide the experi- 
mental verification of the TETRA computer code. Both of these tests involved 
the loss of a single fan blade (or an equivalent unbalance) while running at 
high speed. One of these tests was conducted on a blade-out test vehicle and 
the other on a full engine. Details of these tests and their results are 
given in Sections 4,7 and 4,8. Analytical predictions for both of these tests 
were made using the TETRA computer code and are discussed in Sections 4.4 and 
4.6, The TETRA model for both of these test vehicles was generated using the 
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VAST computer code; this is the primary tool within the General Electric Air- 
craft Engine Business Group for predicting engine system dynamic character- 
istics- This program is an extension of the Prohl-Myklestad method and ad- 
dresses the branched load paths found in real engine structures- The vibra- 
tory behavior of multishaft engine systems can be analyzed, and both synchro- 
nous rotor whirl and distributed linear structural damping are considered. 

The program includes bending, shear deformation, and rotary and gyroscopic 
effects. It computes and prints out the system natural frequencies, displace- 
ments, rotations, and loads for each normal mode, together with the potential 
and kinetic energy distributions in the various components. In addition, VAST 
computes the forced frequency response at each critical speed based on modal 
scale factors derived from the real mode shapes, the energy distributions, and 
user-specified structural damping distributions. 

The VAST program was specifically developed for axisyrametric structures 
typical of most jet engine components. It assumes planar motion with two de- 
grees of freedom per node, i.e., lateral displacement and bending within the 
plane (the gyroscopic forces on the rotor are based on a circular whirling 
motion) ♦ 


4.3 BLADE-OUT TEST VEHICLE VAST ANALYSIS 


The forced response of the test vehicle was predicted using the VAST pro- 
gram. This reponse will provide the basis for assessing the validity of the 
TETRA idealization. 

A schematic diagram of the VAST model of the blade-out test vehicle is 
given in Figure 43. The test vehicle is made up of engine hardware (fan rotor, 
fan frame bearings and fan cases) and supporting slave equipment (drive sys- 
tem, soft mount system, and support structure). Those portions of the model 
representing real engine hardware were taken directly from the VAST model of 
the full engine. New model elements for special slave hardware were generated 
as necessary. In Figure 43, elements identified by circles are spans. A span 
is a structural element made up of small, massless, structural intervals and 
point masses. The structural intervals can represent beams, cones, or springs. 
The point masses can include weight, transverse inertia, and polar moment of 
inertia. The elements identified by triangles are massless springs represent- 
ing bearings, frames, mounts, and supports. The VAST program computes stepping 
or transfer matrices for each span interval and then performs a progressive 
multiplication across each span to obtain equilibrium and compatibility equa- 
tions that are then combined with the system boundary equations to form the 
dynamic equations for the entire system. The blade-out test vehicle VAST model 
consisted of 12 spans, 8 springs (including bearings, frames, and supports), 
and a total of 194 degrees of freedom. The test vehicle was analyzed in the 
vertical direction only. The structure was axisymmetric, as were the primary 
mounts. The secondary supports (spring hangers) were very soft dynamically in 
the vertical direction and provided no support in the horizontal direction. 
These spring hangers were the only major difference between the vertical and 
horizontal planes. Since there was no significant potential energy in these 
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springs in the flexible mode of interest, the critic, .s in the vertical and 
horizontal pla »s were taken to be identical. A mor detailed description of 
the blade-out 2 st vehicle is given in Section 4.7. 
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The two subsystems were then joined in the TETRA system by use of the 
connecting Elements 1 and 2 . A schematic model for the complete TETRA test 
vehicle is shown in Figure 46. This model represents the input used in TETRA 
for the moderate fan-unbalance (linear) analysis. 

TETRA was conducted with 13 casing modes and four rotor modes. The maxi- 
mum frequency used for the casing modes was 75,867 rpm. Fewer than 13 modes 
were initially used to formulate the problem in the TETRA program, but the re- 
sults yielded poor frequency and amplitude matchup for the system resonance 
(4043 rpm) when compared to the VAST analysis for the steady-state forced re- 
sponse case. The 13-mode model gave excellent frequency agreement with VAST, 
but the amplitudes were slightly off. To obtain a better matchup of the No. 

1 bearing load, the damping was adjusted slightly; the casing subsystem Q fac- 
tor was changed from 15.0 to 12.5. This small adjustment was required to com- 
pensate for possible modal truncation errors and to control cost factors in- 
volved with analyzing the test vehicle to higher frequencies. The four rotor 
modes were chosen to give a maximum rotor modal frequency (78,001 rpm) approxi- 
mately equal to that of the highest casing mode used (75,867 rpm). The rotor 
subsystem could have used less than four modes, but no additional cost factor 
was introduced by using four modes since the time-integration step did not 
have to be increased. The casing and rotor subsystem modes were entered into 
TETRA in both the vertical and horizontal directions (similar to VAST, assuming 
a symmetrical system). The connecting elements, the Nos. 1 and 2 bearings, 
were also defined with identical vertical and horizontal spring rates. 

TETRA was initiated at Time » 0 and completed at Time =» 0.32 to give ample 
time for the transient overshoot to be observed. Time steps of 0.00015 sec 
were used for the analysis. This gave 40 time intervals per cycle of vibra- 
tion at the highest frequency mode used in the analysis. The test vehicle had 
only one mass point at which gyroscopic effects were of any importance: the 

fan-rotor stage. This is the only location at which gyroscopic loads were 
introduced in TETRA. 

The time-share TETRA data file for the linear analysis of the test vehicle 
is given in Appendix C, together with a description of key items. 
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The steady-estate response as a function of speed was determined by running 
a succession of solutions initiated at discrete rpm*s to cover the speed range 
of interest. Moderate unbalance levels were used to assure a linear (no rub) 
solution. The resulting, time histories at each speed are given in Appendix D, 
The levels attained after the initial starting transients decayed were used to 
define the steady-state response as a function of speed. The No. 1 bearing 
load response as a function of speed is given in Figure 47. The No, 1 bearing 
housing displacement response as a function of rpm is given in Figure 48, 



Fan Speed, rpm 


Figure 47. Test Vehicle Linear Analysis for Nominal 

Fan Rotor Unbalance - TETRA Results - No. 1 
Bearing Load. 
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Fan Speed, rpm 


Figure 48* Test Vehicle Linear Analysis for Nominal 

Fan Rotor Unbalance - TETRA Results ~ No. 1 
Bearing Housing Deflection. 


The TETRA input for blade-out analysis is similar to that for the moder- 
ate fan-unbalance case. The main difference is the addition of the rub element 
(connecting Element 3, Figure 46) between the fan rotor and fan case to simu- 
late the nonlinear casing constraint on the rotor. This rub element is defined 
in List 2 in the following manner: 

$LIST3 
ITYPE=»3 , 

ILEM-3 , 

JT=i, 35, 

SK=*1000000, 

DBAND»350.0, 

CC=»0, 

$ 

ITYPE»3 identifies the connecting element as being a rub element, ILEM*3 being 
connecting element No. 3 joining points 1 (rotor) and 35 (casing). SK=*1000000 
is the casing ovalization stiffness (Ib/in.) with a dead band (radial clearance 
before hard rub) of 0.350 inch. 

The blade-out analysis was conducted in such a manner as to simulate the 
test vehicle actual speed and decel rate. The restart data file was used from 
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Che linear analysis at 4185 rpm, the speed chac Che blade was released during 
Che cesc. The accual decel race was also simulaced. Figures 49 and 50 show 
Che No. 1 bearing load and deflecCion responses as funccions of cime; Che 
speed crace is also shown. 

The blade-ouC analysis was also conducced ac a sceady scace of 4200 rpm, 
wich and wichouC Che rub elemenC, to demonstrate the difference in transient 
overshoot. Figure 51, wich the rub element activated in the analysis, shows 
the overshoot clearing after 1/2 cycle. With the rub element eliminated from 
Che input data file (Figure 52), a sceady buildup in response occurs and 
reaches a steady-state response after 10 cycles. It is interesting to note 
chat Che rub tends to significantly reduce Che amplitude of the response. 


4.5 BLADE -OOT F.NGTNE VAST ANALYSIS 


As with the test vehicle, a similar VAST analysis was conducted for Che 
complete curbofan engine to evaluate Che TETRA idealization. 

The blade-out test was conducted on a Curbofan engine mounted in a simu- 
laced typical aircraft installation (i.e., inlet, cowling, reverser, and 
pylon). The engine was a two-spool machine wich a high pressure rotor system 
and a low pressure rotor system. For further discussion of the blade-out 
engine test see Section 4.8. The blade-out test engine was analyzed using the 
VAST computer code to predict Che system critical speeds and also the subsystem 
modal data for input into TETRA. The VAST schematic diagram is shown in Figure 
53. The VAST model consisted of 26 spans, 15 springs, and a total of 732 de- 
grees of freedom. 'Five resonances were predicted to be in or near the opera- 
ting range of Che engine (the blade was released at 3366 rpm). They were 
rigid-body modes at 318 and 737 rpm, a pylon-bending mode at 1506 rpm, and 
system modes at 2923 and 3761 rpm. 


4.6 TETRA ANALYSIS OF THE BLADE-OUT ENGINE 


For TETRA) Che engine is represented to be two subsystems, Che rotor and 
casing, joined by four connecting elements: Che Nos. 1, 2, 6, and 7 bearings. 

Subsystem 1, Che casing subsystem, includes Che engine casing, high pressure 
rotor, and aircraft pylon components. The subsystem was analyzed with the 
VAST program to predict Che natural frequencies and mode shapes of the struc- 
ture. The VAST model was supported on its own pylon mounts for the analysis. 
Figure 54 shows a schematic diagram representing this subsystem. A total of 
26 resonant frequencies and modes of the casing subsystem were utilized in 
TETRA; Che highest frequency was 24,268 rpm. 

Subsystem 2, Che low pressure rotor system, was also analyzed wich Che 
VAST program to determine Che system natural frequencies and mode shapes. 

This subsystem was treated as a static (nonrotating) component mounted on two 
very soft springs (100 Ib/in.) without gyroscopic effects. The soft springs 
were used to approximate a free-free system in VAST. A schematic model of 
this subsystem -is shown on Figure 55. Six modes, including two rigid-body and 
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Rotor Speed , rpin Normalized Load 




Time , seconds 

Test Vehicle Blade-Out - TETRA Analysis No. 1 
Bearing Load. 


Figure 49, 










Figure 51. TOTRA Test Vehicle Blade-Out Analysis - Rub Element Included - No. 1 Bearing Loads 
at 4200 rpm. 
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Inlet/Fan Case 



Figure 55. Engine Rotor VAST Schematic, 


the first four bending modes, were used in TETRA with the highest frequency 
being 11,081 rpm. The modal frequencies used in the TETRA analysis for the 
Subsystems 1 and 2 are summarized below: 



1 

353 

1 

55 

2 

731 

2 

73 

3 

1,703 

3 

972 

4 

2,978 

4 

2,421 

5 

5,110 

5 

5,627 

6 

5,370 

6 

11,081 

7 

5,655 



8 

6,485 



9 

6,815 



10 

7,076 



11 

7,852 . 



12 

8,982 



13 

9,644 



14 

10,379 



15 

11,632 



16 

12,665 



17 

13,354 



18 

15,130 



19 

15.469 



20 

18 , 1/2 



21 

18,845 



22 

19,480 



23 

21,616 



24 

22,085 



25 

23,181 



26 

24,268 




The normalized mode shapes for the engine casing and rotor subsystem are given 
in Appendix E. 

The two subsystems were then joined in TETKA by cormj»f: L i ng I 

through 4. The schematic model for the complete TETKA #-ngine is shown in Fig- 
ure 56. This model was used in TETRA for the moderate f/jn«unbal ance analysis. 

As previously noted, TETRA was conducted with 26 casing modes and 6 rotor 
modes. The maximum frequency used in the casing subsystem was 24,268 rpm. 
Attempts using fewer casing modes in TETRA resulted in poor frequency and 
amplitude agreement at the peak response predicted by VAST at 3761 rprni for 
linear forced vibration response. The six rotor modes utilized, with a maxi- 
mum frequency of 11,081 rpm, were deemed sufficient to define the subsystem. 
The current model of 26 cas ing modes and six rotor modes yields excel lent fre- 
quency agreement with the VAST predictions. But, like the test vehicle, the 
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amplitude on the No. 1 bearing loads with Q =* 15.0 gave higher amplitude 
for TETRA compared to the VAST analysis. Similar to the case with the test 
vehicle, a Q of 12.5 in TETRA. would give excellent amplitude agreement on the 
No. 1 bearing loads between the two analyses. Again, as was the case on the 
test vehicle, the change in the Q value is justified to account for possible 
modal truncation effects and the cost factors involved with analyzing the en- 
gine to higher frequency. 

The casing and rotor subsystem modes were entered into TETRA both in the 
vertical and the horizontal directions. Likewise, the connecting elements 
were defined with identical vertical and horizontal spring races. 

TETRA was initiated at Time ® 0 and terminated at Time =» 0.36 to give 
ample time for the transient overshoot to be observed. The time step was 
defined as 50 Usee. This is a conservative time step and is generated by 
giving 40 time intervals per cycle of vibration at the highest subsystem fre- 
quency mode considered for the analysis. Gyroscopic effects were introduced 
at five locations: Joints 1, 7, 8, 12, and 13. 

The TETRA data file for the linear analysis of the engine is given in 
Appendix F. List 1 calls in the casing and rotor subsystems from previously 
created (by VAST) data files. List 2 defines the four connecting elements 
which tie the two subsystems together. And finally. List 3 is the TETRA com- 
mand section of the data file. 

The steady-state response as a function of speed was determined by run- 
ning a succession of solutions initiated at discrete rpm's to cover Che speed 
range of interest. Moderate unbalance levels were used to assure a linear 
(no rub) solution. The resulting time histories at each speed are given in 
Appendix G. The levels attained after the initial starting transients decayed 
were used to define the steady-state response as a function of speed. The 
No. 1 bearing load response as a function of speed is given in Figure 57. The 
No. 1 bearing housing displacement response as a function of rpm is given in 
Figure 58. 

The TETRA input for the fan blade-out analysis on the engine is similar 
to the moderate fan unbalance case with one exception: Che addition of the 

rub element. This rub element is defined in the following List 2: 

200 $ 

201 $LIST2 

202 ITYPE=3 

203 ILEM«5, 

204 JT»1,9, 

205 SK“ 1000000, 

206 DBAND»350,0, 

207 CC=0, 

208 $ 

The engine blade-out analysis was conducted in such a manner to simulate 
the actual sequence of the engine test. The restart data file was used from 
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Figure 57* Engine Linear Analysis for Nominal Fan 
Unbalance - TETRA Results - No, 1 
Bearing Load, 




the linear analysis at 3366 rpm, the speed at which the fan blade was actu- 
ally released. The engine test decel rate was also utilized. Figures 59 and 
60 show the No. 1 bearing load and deflection response as a function of time 
with the decel rate included. 

Figure 59 shows that there was no significant transient load response 
overshoot predicted by TETRA for the blade-out engine test. Figure 60 shows 
that TETRA did predict a transient velocity (or displacement) response over- 
shoot. The magnitude of the predicted transient overshoot was approximately 
45%, The duration of the transient overshoot was very brief and was essen- 
tially over in one cycle. 


4.7 BLADE-OUT VEHICLE TEST DESCRIPTION AND RESULTS 


The blade-out test vehicle consisted of a fan module from a turbofan en- 
gine mounted on soft supports and driven by an electric motor through flexible 
couplings. The fan module consisted of the fan rotor and shaft, the fan frame, 
the fan casing, and the No. 1 and No. 2 bearing supports and bearings. The 
fan module was mounted on a heavy facility support plate mounted from soft 
elastomeric shock mounts and spring hangers. The spring hangers were very 
soft dynamically and were used primarily for vehicle alignment. Except for 
the spring hangers, Che vehicle was totally axisymmetric . The front of the 
fan case was attached with a rubber boot to a larger bell jar. The entire 
vehicle and belljar were filled with helium in order to reduce the horsepower 
required to reach the desired rotor speed. 

The instrumentation used during the blade-out vehicle rig test consisted 
of four No, 1 bearing strain gages and No, 1 bearing vertical and horizontal 
accelerometers. Figure 61 shows the. installation of the instrumentation. The 
strain gages were mounted on the bearing outer race in the hoop direction so 
as to read d 3 mamic load due to the rotating force vector. The four gages were 
90* apart. All of the data were recorded on magnetic tape for future data anal- 
ysis, The test was conducted in two steps. The first step was a mechanical 
checkout of the vehicle up to maximum speed, 4185 rpm, including trim balanc- 
ing the fan rotor, and the second step was the actual bladeout test. Figure 

62 shows the No. I bearing load during the mechanical checkout, and Figure 

63 shows the No. 1 bearing displacement response during the same run. The 
results show that the rotor passes through a response peak at approximately 
4000 rpm. The VAST analysis predicts this peak to be at 4043 rpm. The com- 
parison ot Lest ami VAS’l predictetl bearing load given in Figure 64 shows g<*ii- 
erally reasonable agreement, particularly near the resonance. 

The blade-release mechanism for the fan blade-out test was an explosive 
charge installed in a hole drilled through the root of a single blade just 
above the dovetail. The explosive charge was detonated while running at a 
stabilized speed of 4185 rpm. As soon as the blade was released, the electric 
motor was deenergized, and the rotor coasted down. Figure 65 shows the rpm of 
the rotor as a function of time during the coastdown. All significant unbal- 
ance (released blade plus secondary damage) occurred within 1/4 of a revolu- 
tion. Figure 66 shows the raw signals from two of the No. 1 bearing strain 
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Normalized Foi 



Figure 59. TETRA Engine Blade-Out Analysis - No, 1 Bearing Load 






Normalized Velocity 



Figure 60. TETRA Engine Blade-Out Analysis - No. 1 Bearing Velocity. 
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Figure 61. Blade-Out Test Vehicle Instrumentation 
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Figure 64. Blade-Out Test Vehicle - Test Versus Predicted No. 1 Bearing Load Linear 
Pressure for Nominal Unbalance. 
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Figure 65. Blade-Out Test Vehicle Decel Time Hist 









Normalized Load 


gages 180* apart. Since the gages are measuring hoop strain, they only mea- 
sure an appreciable strain when the load vector is in the direction of the 
gage. When the load vector rotates 180*, there is little compressive strain 
on the gage; therefore, the raw strain signal is asymmetric. In order to 
obtain the true bearing load signal, the high-speed digitizing capability of 
the fast Fourier transform analyzer was used to digitize the data from the two 
gages at 153* and 333®, change the strain to load by applying the nonlinear 
calibration curve, reconstituting the full load versus time signal by changing 
the sign of the 333* strain gage signal, and then combining the two tension 
portions of the two signals. The resulting true load versus time response 
signature is given in Figure 67. Figure 68 shows the same data with an ex- 
panded time scale and the first six cycles of response data. As can be seen 
from the data, there was very little, if any, transient overshoot following 
the blade releane. Figure 69 shows the No. I bearing response (velocity) 
following the blade release. Figure 70 shows the same response with an ex- 
panded time scale. 


Blade Release 



Figure 67. Blade-Out Test Vehicle - Reconstituted No. 1 Bearing Housing Load 
Generated from 153° and 330° Strain Gages. 
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Figure 68. Blade-Out Test Vehicle - Reconstituted No. 1 Bearing Housing Load Generated from 153° and 
333° Strain Gages (Expanded). 







Normalized Velocity 


Blade Release 



Figure 70. Blade-Out Test Vehicle - No. 1 Bearing Response (Expanded). 


4,8 BLADE-OUT ENGINE TEST DESCRIPTION AND RESULTS 


The blade-out engine test was conducted on an operating turbofan engine. 
The test was conducted at the General Electric off-site Peebles Test Operation, 
The engine was mounted in an aircraft pylon and included the inlet, cowling, 
and thrust reverser. The intent of the test was to release portions of two 
fan blades, separated circumferentially by 39*, so that the vector sum of the 
unbalance was equal to. the loss of a single blade. The blade-release mecha- 
nisms were two explosive carriers bolted to the blades at approximately mid- 
span, The two explosive charges were wired to the same detonator circuit so 
that the blades were released simultaneously. High-speed motion pictures 
verified that the two blades were released simultaneously. 

The primary instrumentation for the blade-out engine test was the accel- 
erometers and strain gages on the No. 1 bearing housing. This instrumentation 
was similar to that used during the blade-out rig test (Figure 61), All of 
the vibration and strain signals were recorded on magnetic tape for posttest 
analysis , 
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The blades were released at an engine speed of 3366 rpra. Figure 71 shows 
engine rpm as a function of time following the blade release. Figures 72 and 
73 show the raw strain gage signals following the blade release. These signals 
were reconstituted in the same manner as the rig test bearing strain to give 
the true bearing load signature (Figure 74), As with the vehicle test, tran-» 
sient overshoot after the blade release is small, although a little larger on 
the engine. The No, 1 bearing velocity-response signature is given in Figures 
75 and 76, The displacement response shows a possible small transient over- 
shoot followed by a gradual decrease in response as rpm drops. 


4,9 COMPARISON OF TETRA RESULTS^ FOR THE BLADE -OUT VEHICLE WITH LINEAR 
ANALYSIS AND TEST 


The blade-out test vehicle was analyzed using TETRA for two different 
conditions, i.e,, linear with moderate unbalance and nonlinear with high un- 
balance (loss of a fan blade). The linear analysis results will be compared 
with the VAST linear predictions, and the nonlinear results will be compared 
to the blade-out test results. 

The linear VAST analysis conducted on the test vehicle was summarized in 
Section 4.3. From this analysis, plots of the No. 1 bearing loads and deflec- 
tions as a function of speed were generated for a nominal fan unbalance. The 
TETRA results for this same unbalance were summarized in Section 4.4, This 
analysis considered six different steady-state speed conditions compared with 
the VAST results. As noted in Section 4,4, the Q factor for the TETRA results 
was adjusted slightly from 15.0 (the VAST Q) to 12,5 in order to yield similar 
No. 1 bearing loads. 

The results of this comparison of No. 1 bearing normalized loads and de- 
flections are shown in Figures 77 and 78, respectively. The No. 1 bearing 
loads, although slightly modified to account for modal truncation, show ex- 
cellent agreement for both frequency and amplitude of the response. The No. 

1 bearing deflections show good frequency agreement and a 25% lower response 
for TETRA. Both of these comparisons are made for points on the casing sub- 
system. Figure 79 shows a comparison for the rotor deflections at the fan- 
stage centerline. Again, the frequency response between the two analyses are 
similar; the TETRA amplitude is slightly lower than the VAST predictions (7%). 

Figure 80 is a comparison of the No, 1 bearing transient load signature, 
following the blade release, for both the predicted response (i.e., TETRA 
results) and the experimental data. The TETRA analysis predicts a transient 
overshoot followed by a near-constant-amplitude response. The test data do 
not show a No. 1 bearing transient load overshoot but instead show a gradual 
increase in response over several cycles and then a slow decrease as speed 
falls off (Figure 49) for the speed decay time). 

Figure 81 compares the No. 1 bearing normalized deflection transient re- 
sponse signatures following blade release. The TETRA predicted synchronous 
1/rev deflections are approximately the same on the measured response during 
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Figure 71. Blade-Out Engine Teat Decel Time History, 








Blade Release 



Figure 73. Blade-Out Engine Test - No. 1 Bearing Housing Stress at 333°. 
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Figure 74. Blade-Out Engine Test - Reconstituted No. 1 Bearing Load* 
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Figure 75. Blade-Out Engine Test - No. 1 Bearing Response. 
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Figure 76. Blade-Out Engine Test - No. 1 Bearing Response (Expanded). 
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Figure 78. Test Vehicle Linear Analysis for Nominal Fan Rotor Unbalance - 
TETRA Versus VAST - No. 1 Bearing Housing Deflection. 
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Figure 79. Test Vehicle Linear Analysis for Nominal Fan Rotor Unbalance - 
TETRA Versus VAST - Fan Stage Centerline Deflection. 
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Figure 80. Blade-Out Test Vehicle - TETRA Versus Test 

for Blade-Out Condition - No. 1 Bearing Load. 
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Figure 81. Blade-Out Test Vehicle - TETRA Versus Test for 
Blade-Out Condition - No. 1 Bearing Response. 
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Che first cycle following the blade release. The test data show that the 
response remained approximately constant during the first 30 cycles of vibra** 
tion; however, as the rotor speed decreased the TETRA predicted response shows 
a drastic reduction in amplitude, down to less than 25% of the initial response 
in a few cycles of vibration. Note that the low-frequency content predicted by 
TETRA does not show up in the test data since the accelerometer charge ampli- 
fiers filter out signals below 15 Hz. Figures 82 and 83 compare Che measured 
and predicted No. 1 bearing loads and deflections, respectively, as a fxmction 
of rpm. They show the same story as the time histories: generally reasonable 

correlation between TETRA and the analysis for the loads and very poor agree- 
ment between TETRA predicted displacement response and the measured response. 

It should be noted that the bearing-housing strain gages were calibrated 
in a static-load calibration fixture that accounted for the bearing load dis- 
tribution on the housing. On the other hand, the accelerometer senses response 
at a discrete point and includes out-of-round shell response. It is impossible 
to calibrate the accelerometer to eliminate this effect. TETRA and VAST anal- 
yses assume that plan sections remain round and plane; therefore, it is reason- 
able Co consider that the difference between the measured and predicted perfor- 
mance is due to the method of measurement and not errors in Che prediction 
technique. 


4.10 COMPARISON OF TETRA RESULTS FOR BLADE-OUT ENGINE WITH LINEAR ANALYSIS 
AND TEST COMPARISON 


Comparison of TETRA linear analysis of the engine with the linear VAST 
predictions shows essentially the same results as the comparison for the 
blade-out test vehicle. Figure 84 compares the normalized bearing load for 
Che steady-state solutions versus rpm, and Figure 85 compares the normalized 
deflections. The TETRA and VAST analyses both were conducted using a system 
quality factor Q of 15. If Che TETRA Q had been adjusted to 12.5, like the 
test vehicle data presented in Section 4.9, then the loads would have agreed 
very well, and the TETRA predicted deflections would have been slightly less 
than the VAST predictions. 

Figure 86 compares the normalized bearing load transient response follow- 
ing blade release for the TETRA predictions and the blade-out engine test. 

The signatures look very similar; the maximum load occurs within the first few 
cycles, and the load drops off slowly thereafter as the rpm decreases (see Fig- 
ure 71 for the decel rate). Even chough the signatures look similar, with only 
a small amount of transient overshoot, the TETRA predicted loads are approxi- 
mately 2-1/2 times as high as the measured data. 

Figure 87 compares the TETRA predicted normalized velocity response with 
the test data. The test data show a more gradual buildup to the peak value, 
requiring approximately three cycles of vibration, while the TETRA predicted 
response is approximately twice Che maximum measured values; however, once the 
initial predicted transient overshoot is dissipated, the response levels are 
very similar. Figure 88 compares the normalized loads as a function of rpm 
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Figure 85. Engine Linear Analysis for Nominal Fan Unbalance - TETRA Versus VAST 
No. 1 Bearing Deflection. 



Normalized Load 



Figure 86* Blade-Out Engine - TETRA 
Versus Test for Blade-Out 
Condition - No» 1 Bearing 
Load. 


Figure 87. Blade-Out Engine - TETRA 
Versus Test for Blade-Out 
Condition - No. 1 Bearing 
Response « 
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Figure 88. Blade-Out Engine Test Normalized Load 
Versus rpm. 
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during the decei following the blade release. Figure 89 shows similar curves 
for Che normalized deflections. Figures 88 and 89 both show that Che response, 
a function of rpm, is very similar between TETRA predictions and measured data 
except for the differences in magnitude and transient overshoot previously 
noted. 


4.11 DISCUSSION OF TETRA VERIFICATION RESULTS 


The TETRA linear analyses agreed reasonably well with the VAST linear 
analyses, especially in terms of response frequencies, for both the test vehi~ 
cle and the engine. Both analyses were dependent on the number of modes used 
to define the subsystems; the engine analysis were the more critical due to 
the increased complexity of the model. The discrepancies in the TETRA results 
are in the direction Chat might be expected due to modal truncation; i.e., fre- 
quency and loads are slightly too high, and deflections are slightly too low. 
Table IX is a comparison of the complexity of the vehicle and engine models. 

The comparison of the nonlinear results showed a much larger disparity 
between the test vehicle and the engine. The vehicle showed reasonable agree- 
ment on the bearing loads, but the predicted deflections were much lower than 
test. For the engine, predicted loads were more than twice the measured loads, 
but the deflections showed reasonable agreement following the initial tran- 
sient overshoot. For both the vehicle and the engine, TETRA predicted more 
transient overshoot than was observed during the blade-out tests. 

It was apparent from preliminary analyses and the correlation studies 
that the number of modes used to represent the subsystems is a key to the 
accuracy of the TETRA response. No proven guideline is currently available 
to define what a reasonable number of modes should be; however. Table IX may 
give some insight. Table IX shows that, even though the number of casing 
modes used in the engine analysis was twice that used in the vehicle analysis, 
the ratio of the highest casing frequency to the maximum operating frequency 
was only one-third of that for the vehicle. Likewise, the ratio of the number 
of degrees of freedom to the number of modes used was twice as high for the 
engine although it is not clear that the number of subsystem degrees of free- 
dom is a relevant parameter. It was apparent from the data that the test vehi- 
cle analysis more closely represented both the linear VAST analysis and the 
test results; this would imply that improvements in the engine predictions 
could have been obtained by using similar ratios of maximum subsystem fre- 
quency to operating speed and/or number of degrees of freedom to the number 
of modes used. Incorporating these criteria into the engine analysis would 
have increased the computational cost to perhaps prohibitive values: approxi- 

mately $3500 for 0.5 seconds of data; with the current analyses, each cost 
approximately $600. This may not be a prohibitive value for a design analysis, 
if there is confidence in the results and it would not require repetition to 
evaluate the effect of other combinations of modes on the results. A possible 
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Figure 89. Blade-Out Engine Test Normalized Deflection. 
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Table IX, TETRA Model Complexity. 


Parameter 

Vehic le 

Engine 

Max. Operating Speed 

4,185 rpm 

3,366 rpm 

Max , Cas ing Frequency 

75,867 rpm 

24,268 rpm 

Max. Rotor Mode 

78,001 rpm 

11,081 rpm 

Max, Casing Freq. 
^ , Oper . Speed 

18.1 

7.2 

. Max Rotor Freq. 

Max. Oper. Speed 

18.6 

3.3 

No. DOF, Casing Subsystem 

170 

590 

No. DOF, Rotor Subsystem 

38 

190 

No. of Subsystems 

2 

2 

_ Casing DOF 

Ratio — -a-; — 

No. Modes 

13.1 

22.7 

„ , . Rotor DOF 

No. Modes 

9.5 

31.7 

No. of Connecting Elements 

2 

4 

Cost to Run 0.5 second 

$600 

$630 


way of increasing the accuracy would be to break the engine down into more and 
simpler subsystems such as the low pressure rotor, the high pressure rotor, the 
high pressure casing, the low pressure casing, and the pylon. Each of these 
systems would require fewer modes to get a high ratio of subsystem frequency to 
operating frequency and also to achieve a low ratio of degrees of freedom to 
number of modes used. The integration time step would have to be small to 
account for the higher frequency, and the analysis would be more complicated 
with more subsystems and more connecting elements, so the true effect on cost 
is unknown. 

Other factors which may have had an effect on the comparison of predicted 
response and test results are: 

« System damping at high amplitudes 

e Tip clearance and casing ovalization stiffness 
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m Accuracy of measured loads and deflecCions 

m System complexity. 

It has been observed in other studies that the Q factor has a small 
effect on the nonlinear response. The tip clearance and ovalization stiff- 
ness were determined by measurement and finite-element analysis of the casing 
structure. In the past, small variations of these parameters have been seen 
to have a minor effect on the results. As was mentioned in Section 4,9, there 
are some inherent inaccuracies in the deflection measurements which could ac- 
count for some of the differences between the predicted and measured response 
of the blade-out test vehicle. 

From these comparisons between. TETRA calculated responses and the test 
results from both test vehicle and engine blade-out tests, the following 
observations are made: 

- e TETRA shows reasonable agreement with both linear analysis and non- 
linear test results when a sufficient number of modes are included. 

• TETRA tends to predict more transient overshoot than was observed 
during both blade-out tests but significantly lower overshoot than 
predicted by linear analysis (no rub). 

e Transient overshoot is insignificant in the test data, 

# Modal truncation can have a significant effect on the accuracy of 
TETRA computed response. 

e Other . real-world factors, difficult to quantify, affect the response 
of real engines. Some of the more important ones are plastic de- 
formation and the effect of high impulsive loading on structure 
damping and stiffness. 
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5.0 CONCLUSIONS AND RECOMMENDATIONS 


A comparison of analytical predictions with test data for the two 
configurations tested was given in Sections 4.9 and 4.10. The essence of 
these comparisons and other observations is: 

• Loads were more accurately measured in test than were displacements 
and, therefore, are taken as the prime reference for comparison. 

• Analysis and test agree well if a sufficient number of subelement 
modes are used. 

• The test vehicle analytical results agree very well with test, 
but the engine results compare poorly. The latter used sub- 
element modes up to a maximum frequency of roughly one-third 
that used for the test vehicle. 

• Discrepancies between analysis and test depend strongly on the 
maximum frequency of the subelement modes relative to the maximum 
operating frequency. 

• Test results show relatively little transient overshoot in loads; 
the analysis shows about 30%. 

To attempt a rough quantification of results, a comparison is given below 
between predicted loads and test results both for the transient overshoot and 
the steady-state values. Some license was taken in developing this table since 
the determination of transient and steady-state values is subjective, espe- 
cially since no steady-state condition was truly obtained. While it is readily 
conceded that other parameters (such as the number of subelements, etc.) may 
have as important an influence, the table of comparisons is made only on the 
basis of the ratio of the maximum subelement modal frequency to the maximum 
rotational frequency of the test. 


Max Modal Freq. No . of 
Configuration Rotational Freq. Modes 


% Difference 
Steady-State 
Load 


Test-Analysis 

Test 

» Transient 
Load , % 


X 100 


Test Vehicle 

18 

D 

8% 

-33% 

Engine 

7 

26 

144% 

131% 
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The test vehicle results show very good correlation, especially for 
steady-state response, but Che engine shows poor. Also, the analysis tends to 
predict greater overshoot than actually occurs. An implication is that even 
better results from the test vehicle might be obtained, especially for the 
transient response. It is not clear chat the only factor producing the 
difference noted between analysis and test is related solely to the maximum 
subelement modal frequency. As mentioned previously, how the system is 
modeled may be equally important. For instance, the engine was modeled as two 
subelements, one representing the LP rotor and the ocher representing a 
combination of basically soft and stiff elements such as the pylon, casings, 
frames, and HP rotor. Perhaps better results are possible if the pylon, 
casing, and HP rotor are each represented as a subelement. 

Also, it may be possible to obtain better results if a different sub- 
element modal representation is used. In the approach used in TETRA, free- 
free modes are used to define the subelement (except for elements tied to 
ground, for which grounded boundary conditions are used). It seems reasonable 
to believe that faster convergence and fewer modes would be required if the 
subelement modes were taken to be pinned-pinned at Che connection points, or 
tied to ground through springs having the support structure spring rates. 

These modes could be closer in frequency and shape to the dominant modes of 
the combined system and, therefore, could require fewer modes to obtain good 
results. Another complicating factor is the influence of fan tip rub which 
produces a secondary load path and essentially a "stiff er" system. Hence, 
a greater number of modes would be required for this condition than one would 
expect to use to obtain good results for the linear case. 

It is evident from the results of these discussions that criteria and 
guidance are needed if the analyst is to precede with assurance that he will 
predict accurate loads and stresses in a cost-effective manner. These guide- 
lines and criteria are not at hand and must be developed as part of this 
approach. The question of minimizing truncation errors for complex systems 
involving nonlinear stiffening elements is difficult.. Furthermore, the most 
effective modeling approach must be identified. 

Based on these discussions, the following list of areas requiring further 
study was compiled: 

1. Truncation Errors - Need to develop criteria for determining the 
number of subelement modes required to produce acceptable results, 
including the effect of nonlinear elements. 

2. Modeling Guideline - Need to develop criteria for guiding the 
analyst in the selection of the best subeleraent breakdown. Criteria 
should include a balance between increasing number of subelements, 
diminishing time increments, and improved accuracy consistent with 
cost. 
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3. Component Modes ^ The current analysis uses free-free subelement 
modes at discrete points, usually at bearings or frames. The use of 
subeleraent normal modes pinned-pinned at these points, or sprung by 
elements having stiffnesses equal to the connecting element stiff- 
ness, may produce better results and require fewer modes. These 
approaches should be explored. 

4. Transient Overshoot - The criteria and modeling techniques studies 
discussed above, must also address the problem of predicting over- 
shoot or Lack thereof. This may also include criteria defining 
the required time increments. 

5. Cost - The cost of running complex, real-engine systems using the 
large number of modes indicated by the results and the large number 
of time increments required to obtain good results can be prohibi- 
tive. We need to study methods for reducing computer time and cost. 

6. Instrumentation - Accurate measurements of displacements and loads 
under high-amplitude and duration testing is difficult, especially 
when local shell deformations occur. A further study of this area 
is needed. 

7. Experience - More experience in the use of this approach is required 
on a variety of complex problems to explore its full capabilities. 

The work and methods developed under this contract provide a useful tool 
that was demonstrated to give good prediction of blade-out test results. 

This study also showed that the use of this tool is not straightforward, and 
the analyst must be provided with criteria and guidelines to assure good re- 
sults. Developing these guides is the next logical step for making this 
method a useful design tool. 
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APPENDIX A 


BLADE-OUT VEHICLE SYSTEM NORMAL MODE SHAPES 


The mode shapes for the first three resonances predicted by analysis of 
the blade-out test vehicle are given in Figures 90, 91, and 92, The first 
two are rigid-body modes, and the other is a system flexural mode* 


Preceding page blank 
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Normalized Deflection 



Figure 90. 


BladeHDut Test Vehicle Predicted Mode Shape at 
387 rpm. 
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Normalized Deflection 



Figure 91 . Blade-Out Test Vehicle Predicted Mode Shape at 
981 rpm. 
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Normalized Deflection 



Axial Length 


iT 

Figure 92. Blade-Out Test Vehicle Predicted Mode Shape at 
4043 rpm. 
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APPENDIX B 


BLADE-OUT TEST VEHICLE SUBSYSTEM NORMAL MODE SHAPES 


Figures 93 through 105 show the natural frequencies and normalized mode 
shapes predicted for the casing subsystem. Figures 106 through 109 present 
the frequencies and normalized shapes for the rotor subsystems, 
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Figure 100. Blade-Out Test Vehicle, Subsystem 1 (Casing), 36, 388 rpm 
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Figure 101« Ulade-Out Test Vehicle, Subsystem 1 (Casing), 39,995 rpm 




Figure 102, Blade-Out Test Vehicle, Subsystem 1 (Casing), 41,284 rpm 
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Figure 103. Blade-Out Test Vehicle, Subsystem 1 (Casing) , 53,107 rpm. 















69,016 rpm. 
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Figure 105, Blade-Out Test Vehicle, Subsystem 1 (Casing), 75,867 rpm. 
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Figure 107. Blade-Out Test Vehicle, Subsystem 2 (Rotor), 146 rpm 
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Figure 108* Blade-Out Test Vehicle, Subsystem 2 (Rotor), 17,626 rpm 






APPENDIX C 


TETRA DATA FILE FOR THE BLADE-OUT VEHICLE LINEAR ANALYSIS 


Figure 110 is a listing of the time-share TETRA data file for the linear 
analysis of the test vehicle. 

LISTI (Lines 100-760) - Lines 200 and 201 are control statements, for the 
casing sybsystem, that call in the model data in the vertical and horizontal 
directions respectively. The model files were previously created by the VAST 
computer code. For Figure 110, the same model data is used in either direc- 
tion. The number 15.0 is the Q of the casing subsystems. Lines 220 and 224 
are the calling statements for the rotor subsystem data files. 

LIST2 (Lines 770-1000) - These lines define the connecting elements uti- 
lized to join the casing and rotor subsystems. 

Lists 1 and 2 define the boundary conditions of the test vehicle. List 3 
is the command section of the data file. The key items on this list are the 
following: 

Line No. 

1090 Definition of the beginning speed (4100). 

1120 Magnitude of the unbalance (gram-inches) time at 

initiation of unbalance * 0.00004 seconds with the 
unbalance located at Joint I . 

1130 Definition of joints where gyroscopic effects need 

to be considered. 

1150-1330 Definition of joints and directions requested for 

the plotting output file. 
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\iS $LIST1 

113 KAKE=63HBACH-GALLARB0 DATE 8*83 CHARGE 18513 
123 ADDRE3-83H 533 BSHT KAIL DROP K7l EXT 4279 
132 IDENT1--63HKHIRLIGIG TEST VEHICLE TETRA SETUP STAHB HCU3IHG 
140 IB£KT2=80KCHECK CASE VERT % HOR2 KODEL -GROUNED CASIHG 


200 S3arl)=7.0i0.0!l5.0»15r 
210 3i0t0i0*15.0il5> 

220 If0f0f0>0il8t 

225 2>0*0>0>0>18 
760 $ 

770 $LIST2 
780 ITTPE--5. 

790 ILEK=lr 
800 UT:4>2i 
305 rS=1.183£6. 

810 Z3M.183E6> 

850 IDAKPslr 
360 8ELEK--15. 

870 8FRES=67.3. 

330 $ 

890 $LIST2 
900 mPEsS. 

910 ILEK--2. 

920 JT=20.3» 

925 TS=3.4E6. 

930 23s3.4E6t 
970 IDAKP--l» 

980 8EL£H*15f 
990 8FREG=67.3f 
1000 ( 

1010 6LIST3 
1020 ISTARTi0. 

1030 DELTA=0.000015> 

1040 TFIKALs0.32i 
1050 IPRKUL=500f 
1060 IPLKULS10. 

1070 IRCTI-1» 

1080 BEGTIK--0, 

1090 BEGRPKS4100. 

1100 TRHI3(lfl)= 

1110 .32.0. 


1120 UHBAL(1.1)=. 00004.1. .0. 

1130 GTROd.Ds 
1140 1.3.28E5. 

1150 NPBll.Ds 
1160 1.1. 

1170 35.1. 

1180 2>lf 
1190 4.1,- 
1200 3.1. 

1280 KEPD(1,1)> 

1290 1,4,1, 

1295 1,4,3. 

1310 2,20.1, 

1315 2,20,3, 

1330 $ 


Figure 110. Listing of TETRA Input Data File. 
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APPENDIX D 


BLADE--OUT VEHICLE TETRA RESULTS FOR MODERATE UNBALANCE 


Steady-state response as a function of speed was determined by running a 
succession of solutions initiated at discrete rpm's to cover the speed range 
of interest. Moderate unbalance levels were used to assure a linear (no rub) 
solution. The resulting time histories at each speed are given in this 
Appendix. The levels attained after the initial starting transients decayed 
were used to define the steady-state response as a function of speed. 

The individual TETRA time histories with moderate unbalance are given 
in Figures 111 through 122. These plots show No. 1 bearing loads and deflec- 
tions (vertical direction) at each steady-state speed point as a function of 
time. In each case, the unbalance was initiated at Time » 0.0004 second. 
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Figure 112. No, 1 Bearing Deflection at 3500 rpm. 




Figure 113. No. 1 Bearing Load at 3800 rpm. 








184 


Figure 116. No. 1 Bearing Load at 4000 rpm. 
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Lng Load at 4050 rpm. 
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Figure 121* No. 1 Bearing Load at 4400 rpni. 
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Figure 122* No* 1 Bearing Deflection 4400 rpm. 


APPENDIX E 


BLADE-OUT ENGINE NORMALIZED SUBSYSTEM MODE SHAPES 


The nonaalied mode shapes for Che casing subsystem are presented In 
Figures 123 through 148. Normalized mode shapes for Che rotor subsystem follow 
in Figures 149 through 154. 
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Axial Length 

Figure 123. Subsystem 1 (Casing) Normalized Mode Shapes at 353 rpm. 
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Axial Length 

Figure 125. Subsystem 1 (Casing) Normalized Mode Shapes at 1703 
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Figure X26. Subsystem 1 (Casing) Normalized Mode Shapes 
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Figu 
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Axial Length 

Figure 138. Subsystem 1 (Casing) Normalized Mode Shapes at 12,665 rpm. 
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Figure 140. Subsystem 1 (Casing) Normalized Mode Shapes at 15,130 rpra. 
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Figure 141* SubBystem 1 (Casing) Normalized Mode Shapes at 15,469 
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Figure 145. Subsystem 1 (Casing) Normalized Mode Shapes at 21,616 rpm 
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APPENDIX F 


TETRA DATA FILE FOR THE BLADE-OUT ENGINE LINEAR ANALYSIS 


Figure 155 is a listing of the TETRA data file for the engine. List 1 
calls in Che casing and rotor subsystems from previously created (VAST) data 
files. List 2 defines the four connecting elements, and List 3 is the command 
section. 
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4*>/t ±t 

kLO 9L10U 

130 KAK£=60KBACH-SA'lLARBO DATE 3-30 CHARGE 13512 
142 ADBRE3S62H 522 B3RT RAIL DROP K7l EXT MTi 
152 IBEKTl=fc2H CF6-52 «I1C OiGIHE TETRA CHECKOUT 
132 IB3{T2=42KCKECK CASE VERT AND HORIZ RODEL OKLT 


172 SSdilJs 

182 7i2>2i2tlS.2il5> 

185 8r2i2i2)15a2tl5i 

132 I>2i2i2t2>18* 

135 Z>2i2i2i2fl8i 

222 i 

212 $LI3T2 

222 ITYPE=5. 

232 ILER=1» 

111 

245 Y3s4.2£3i 
252 ZS-4.2E3I 
232 IBARP^lt 
272 QELER=15.2» 

282 QFRES=58.5f 
232 $ 

A/»/T #1 f/^TO 

312 ITYPE=5» 

322 ILER=2. 

332 JT=3»11» 

335 YS-3.331E3! 

342 23=3.331E3i 
350 IBAKP=1> 

332 QELERM5.2. 

372 QFR£a*-58.5. 

382 $ 

332 $LIST2 
422 ITYFE-5* 

412 IL£R-3» 

422 vlT=4il7i 
425 YS=3.115E3i 
432 ZS--3.115E6. 

442 IBARpsh 
452 8£L£R-15,2» 

432 0FRE8:58.Si 
472 $ 


482 $LIST2 
432 ITYPE=5i 
522 ILEH*-4, 

512 JT=5il8» 

515 YS=2.543E3. 
522 ZS=2.543E3i 
532 IDARP^i 
542 &ELB:1S.2> 
552 9FREa-58.5> 
532 i 
572 $LIST3 
582 ISTART=2» 

532 DELTA--.22225» 
322 TFIHAL:.33i 
312 IPRHUL:2S2i 
322 IPLRUL--5» 

332 IRCTMi 
342 BEGTIR=2f 
350 BECRPR=3722f 
332 TRHISdfl): 
372 .33>2» 


OGI/ Uli&MLiifi/ 

332 2.222251 1) 
722 GYRCd.l'/= 
712 1 >3282201 
722 7>32252> 

/w0 

742 12i27322i 
752 13>27222> 
732 NPDdfl)= 
772 Iflf 
732 2»r> 

732 3il> 

322 4il> 

812 5>1> 

822 3>1> 

832 7fli 
842 3>li 

AM ^ A 4 

ojr i^fll 

332 13>1) 

872 12tl> 

882 ll>li 
332 14>lf 
322 15>1> 

312 13>1> 

322 13>1> 

AA.T « 

iov ut^ii 

9 tit 4 

'fV Ufkf 

352 18>1> 

332 23>li 
372 24>li 
382 25>1> 

332 21>1> 

1222 22>li 

1212 3>1> 

1222 i«PBd>ll 
1232 1>2>1> 

4atit ^ A 4 


1232 4|5>1> 
1272 $ 


Figure 155, TETRA Data File Listing for Engine. 
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APPENDIX G 


BLADE-CUT ENGINE TETRA RESULTS FOR NOMINAL UNBALANCE 


The steady-state response as a function of speed was determined by run- 
ning a succession of solutions initiated at discrete rpm*s to cover the speed 
range of interest* Moderate unbalance levels were used to assure a linear 
(no rub) solution* The resulting time histories at each speed are given in 
this Appendix* The levels attained after the initial starting transients 
decayed are used to define the steady-state response as a function of speed* 
The individual TETRA time histories with moderate unbalance are given in 
Figures 156 through 163. These plots show the No. 1 bearing loads and deflec- 
tions (vertical) at each steady-state speed point as a function of time. In 
each case, the unbalance was initiated at Time » 0.00025 second. 
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Figure 156. No. 1 Bearing Load at 2950 rpm 
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Figure 157. No. 1 Bearing Deflection at 2950 rpm. 
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Figure 158, No. 1 Bearing Load at 3600rpm. 
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Figure 159* No. 1 Bearing Deflection at 3600 rpm. 
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Figure 160. No. 1 Bearing Load at 3750 rpm. 
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Figure 161. No. 1 Bearing Deflection at 3750 rpm. 
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Figure X62. No. 1 Bearing Load at 4000 rpm. 
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